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INTRODUCTION 

The  NAVAIRDEVCEN  (Naval  Air  Development  Center)  is  a  large  user  of  vehicle  and  propul¬ 
sion  design  and  performance  computer  codes  and  is  constantly  seeking  to  improve  their  efficiency 
and  flexibility,  A  significant  number  of  these  codes  are  dependent  on  the  use  of  input  tabular  data 
sets.  Quite  frequently  these  data  sets  are  initially  received  in  a  format  incompatible  with  direct  use 
in  these  codes,  resulting  in  a  time  consuming,  error  prone  transformation  task.  To  circumvent  this 
problem,  development  of  a  rapid  data  transformation  code  was  undertaken.  The  impetus  for  this 
effort  was  the  need  to  prepare  for  a  planned  substantial  Increase  in  analyses  of  various  aircraft  and 
propulsion  systems. 

This  present  report  describes  a  code  based  on  the  use  of  an  interactive  graphics  system  that 
permits  direct  creation  of  digital  tabular  data  sets  from  material  in  graph  form,  utilizing  a  Tek- 
tronics  4015  graphics  terminal,  digitizer  tablet  and  hardcopy  unit.  In  addition  the  user  may  edit 
and  correct  these  data  directly  from  the  digitizer  tablet  or  from  the  graphics  display  screen  using 
cursor  cross  hairs  and  tablet  commands.  This  code,  entitled  TIGS  (Table  Plot  Interactive  Graphics 
System)  was  developed  using  the  NAVAIRDEVCEN  CDC  6600/Cyber  175  computer  facilities. 

A  user's  guide  for  this  code  is  shown  in  Appendix  A.  A  Fortran  listing  of  the  TIGS  code  is  shown 
in  Appendix  B. 


DISCUSSION 


CODE  DEVELOPMENT 

The  TIGS  code  was  developed  as  a  general  purpose  computer  tool  to  permit  the  user  to 
prepare  and  edit  tabular  data  sets,  using  interactive  graphics,  prior  to  use  in  other  computer 
codes.  The  tabular  data  sets  may  represent  a  functional  relationship  between  a  dependent  variable 
and  several  independent  variables,  an  example  of  which  is  shown  in  Figure  1.  In  this  figure  FXYZ 
is  the  dependent  variable  and  is  a  function  of  the  independent  variables  X,  Y,  and  Z.  The  basic 
output  of  the  code  is  graphical  plots  on  a  Tektronics  4015  type  of  storage  tube  graphics  terminal 
along  with  a  computer  file  consisting  of  the  digital  tabular  data  representation  of  that  plot.  These 
digital  tabular  data  are  suitable  for  use  in  nearly  all  of  the  vehicle  and  propulsion  design  computer 
codes  used  within  the  Aircraft  and  Crew  Systems  Technology  Directorate  at  the  NAVAIR¬ 
DEVCEN.  Further  details  of  the  tabular  data  output  format  are  discussed  in  the  user's  guide 
Appendix  A  and  in  reference  (a).  While  the  TIGS  code  is  a  stand-alone  interactive  system,  the 
graphical  executive  portion  of  the  code  may  be  used  in  conjunction  with  any  other  user  written 
code.  In  effect  this  flexibility  permits  the  user  to  interactively  prepare  and  edit  data  which  in 
turn  is  passed  to  the  user's  code.  Experience  in  using  TIGS  has  shown  that  the  time  required  to 
prepare  data  for  use  in  the  vehicle  and  propulsion  design  codes  has  been  reduced  by  a  factor  of  10. 

HARDWARE  REQUIREMENTS 

The  TIGS  code  is  specialized  in  that  it  was  written  for  a  CDC  6600/Cyber  175  computer 
system  using  a  1200  baud  line  under  the  CDC  telex  time  sharing  system.  Graphical  implementa¬ 
tions  are  provided  by  a  Tektronics  model  4015  terminal  with  the  enhanced  graphics  option.  A 
large  Tektronics  tablet  may  be  employed  in  the  digitization  process  along  with  a  model  4631  hard¬ 
copy  unit.  The  TIGS  system  could  be  modified  for  use  with  other  graphics  systems.  Figure  2 
shows  a  typical  TIGS  hardcopy  plot. 
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SOFTWARE  OVERVIEW 

The  TIGS  code  is  comprised  of  seven  basic  modules  using  the  standard  utility  Tektronics  re¬ 
lease  3.2  software  compiled  under  Fortran  IV.  TIGS  uses  the  Cyber  segmentation  loader  requiring 
about  40000  octal  memory  locations  to  execute.  The  segmentation  setup  consists  of  seven  modules 
described  below.  The  information  flow  between  these  modules  is  represented  by  Figure  3. 

TIGS  is  the  main  executive  module  that  controls  the  input  and  output  and  interplays  with 
the  graphics  executive. 

TABR  contains  the  code  to  input  and  output  the  digital  data  in  the  required  format. 

TIGPPR  is  the  graphics  executive  module.  This  module  controls  the  graphical  input  and  per¬ 
mits  the  user  to  interact  with  the  graphical  screen  and  digital  tablet  controlling  data  point  values, 
plot  sizes,  curve  options,  titles  and  scaling. 

The  TIGPPR  module  performs  these  functions  through  connections  to  other  segmentation 
modules  GETVAL,  LOPTIM,  LABEL,  and  DRAWIT. 

GETVAL  is  used  to  input  data  points  either  from  the  graphics  screen  or  the  digitizer  tablet. 

LOPTIM  implements  the  axes  scaling  and  grid  options  selected  by  the  user. 

LABEL  uses  the  data  values  to  compute  the  axes  tic  marks  and  other  data  related  to  fitting 
the  plot  on  the  graphical  screen. 

DRAWIT  processes  the  scaling,  axes,  along  with  other  plot  data  and  generates  the  commands 
that  draw  the  vectors  on  the  graphic  screen. 

There  are  two  basic  operating  modes  in  the  graphics  executive:  creation  and  correction.  In 
the  creation  mode  a  digital  data  file  is  created  using  the  cross  hair  cursor  either  directly  from  the 
Tektronics  screen  or  from  the  digitizer  tablet.  Commands  from  the  screen  are  implemented  by 
first  positioning  the  cross  hairs  and  then  keying  a  single  letter  indicating  the  command.  Commands 
from  the  digitizer  tablet  are  implemented  in  two  steps:  first  the  command  code  letter  is  keyed 
using  a  tablet  command  menu;  second,  the  coordinate  position  going  with  the  command  is  keyed 
at  the  desired  position.  From  either  the  screen  or  the  tablet,  the  graphics  executive  receives  the 
command  and  coordinate  position.  The  commands  received  by  the  graphics  executive  are  generally 
used  in  three  different  ways: 

1 )  add,  delete  or  change  a  coordinate  point 

2)  change  a  graphics  executive  switch  from  off  to  on  or  on  to  off 

3)  control  the  size  and  view  of  the  graphical  plot. 

Some  commands  available  on  the  screen  can  not  be  used  on  the  digitizer  tablet.  A  more 
detailed  discussion  of  these  commands  is  found  in  the  user's  guide.  Appendix  A. 

USER  EXPERIENCE 

The  TIGS  system  has  proved  to  be  a  very  powerful,  flexible  tool.  Task  cost  reductions  of  10 
to  1  have  been  shown  to  date  using  TIGS  to  prepare  tabular  inputs  for  other  codes. 
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FIGURE  2.  TIGS  EXAMPLE  PLOT 
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FIGURE  3.  TIGS  SEGMENTATION  MODULES 
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CONCLUSIONS 

An  interactive  graphical  code  system  has  been  developed  which  is  capable  of  rapid  transfor¬ 
mation  of  graphical  information  into  tabular  data  formats  which  are  compatible  with  the  input 
requirements  for  a  large  variety  of  in-house  programs. 
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A.1  INPUT  CONSIDERATIONS 

The  TIGS  code  will  permit  a  direct  creation  of  a  data  file  from  screen  and  or  tablet  commands. 
In  addition,  existing  table  data,  input  as  file  TAPE1,  may  be  edited  and  corrected.  In  either  case  a 
new  table  data  source  with  corrections  is  produced  as  an  output  on  the  TAPE?  file.  The  format  of 
files  TAPE1  and  TAPE?  is  the  same.  In  the  next  section  this  format  is  illustrated. 

A.2  TABLE  DATA  FORMAT 

The  table  data  may  represent  a  dependent  (output)  parameter  as  a  function  of  1,2,  or  3 
independent  (input)  parameters.  The  basic  method  for  inputting  these  tables  is  described  in 
reference  (a).  However,  for  the  sake  of  completeness,  the  user's  guide  portion  of  reference  (a)  has 
been  extracted  and  duplicated  herein  and  includes  those  modifications  introduced  since  its  initial 
publication.  This  information  is  shown  in  Table  A-1. 

A.3  EXAMPLES 

Card  input  data  set-ups  for  three  different  examples  are  illustrated  as  follows: 

Example  1  (Drag  coefficient  as  a  function  of  Mach  number) 

The  dependent  variable  is  drag  coefficient  and  the  independent  variable  is  Mach  number. 

Figure  A-1  illustrates  the  graphical  relationship.  This  is  a  one  parameter  table  look-up  so  the 
other  two  parameters  are  dummies.  Table  A-1  shows  the  card  set-up  for  this  example.  The 
EOT  (end  of  table)  parameter  label  terminates  the  data  for  this  table. 

Example  2  (Drag  coefficient  as  a  function  of  Mach  number  and  lift  coefficient) 

The  dependent  variable  is  drag  coefficient  and  the  independent  variables  are  Mach  number 
and  lift  coefficient,  illustrated  in  Figure  A-2.  This  is  a  two  parameter  table  look-up  so  that 
the  third  parameter  is  a  dummy.  Table  A-2  shows  the  card  set-up.  In  Table  A-2  the  last 
Mach  parameter  data  repeats  the  previous  Mach  parameter  data.  In  this  situation,  the  last 
Mach  parameter  data  card  can  be  omitted.  As  a  general  rule,  whenever  the  data  on  the  X 
parameter  axis  as  shown  in  Figure  A-1,  is  repeated,  then  the  X  parameter  data  card  need 
not  be  repeated. 

Example  3  (Drag  coefficient  as  a  function  of  Mach  number,  lift  coefficient  and  CG  location) 
The  dependent  variable  is  drag  coefficient  and  independent  variables  are  Mach  number,  lift 
coefficient,  and  CG  location,  illustrated  in  Figure  A-3.  Table  A-3  shows  the  card  set-up  for 
this  three  parameter  example.  Note  that  the  input  card  set-up  is  symmetrical  in  that  each  CL 
parameter  data  card  begins  data  for  each  CG  parameter. 

A.4  LIMITATIONS 

The  TIGS  system  as  presently  written  is  limited  to  a  maximum  of  30  curves  per  plot,  150 
points  per  curve,  or  a  total  of  300  points  per  plot.  For  example,  a  plot  with  10  curves  could  be 
described  with  5  curves  using  40  points  per  curve,  and  the  remaining  5  curves  using  20  points  per 
curve. 

A.5  INTERACTIVE  PROMPTING 

The  TIGS  code  has  been  designed  to  prompt  the  user  in  supplying  information  in  the  correct 
format.  Selection  of  the  baud  rates  compatible  with  available  transmission  lines  is  possible.  After 
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TABLE  A-l 

DATA  INPUT  INSTRUCTIONS 


1  Table  re-ference  number  of  table  look-up  function;  table  title  or  descriptive  1X,I47A10 


information 

2  4  character  identifier  (user  selected)  used  to  identify  the  third  independent  A4,I3,3X,7F10.0 

variable.  If  table  look-up  has  2  or  less  independent  variables,  use  a  dummy 
identifier;  the  number  of  values  of  the  third  independent  variable  (must  be 
less  than  100.);  values  of  the  third  independent  variable  arranged  in  ascending 
order. 

2a,b,etc.  Continuation  of  third  independent  variable  array,  if  required  10X,7F10.0 

3  and  All  remaining  cards  have  the  same  format  as  card  2,  2a,  b,  etc.  The  item 

following  which  distinguishes  the  card  types  is  the  value  of  the  independent  variable. 

The  4  character  identifiers  of  each  independent  variable  must  not  be 
identical.  The  4  characters  of  each  independent  variable  card  (after  the  title 
card)  are  user  selected.  The  card  order  of  each  independent  variable  is  sig¬ 
nificant.  The  first  four  cards  with  respective  independent  variables  are  as 


follows; 

Card 

Definition 

2,2a,b,etc. 

third  independent  variable, 
identifier  and  values 

3,3a,b,etc. 

second  independent  variable. 
Identifier  and  values 

4,4a,b,etc. 

first  independent  variable, 
identifier  and  values 

5,5a,b,etc. 

dependent  variable,  Identifier 
and  values 

The  remaining  input  cards  use  these  same  identifier  values  as  input  above.  On 
cards  4,4a,  b,  etc.  and  5,  5a,  b,  etc.  are  the  dependent  and  first  independent 
variable  values  along  the  line  given  by  the  first  value  of  the  second  Independent 
variable  and  in  the  plane  of  the  first  value  of  the  third  independent  variable. 

Cards  with  the  same  respective  identifier  value  cards  4,  4a,  b,  etc.  and  5,  5a,  b, 
etc.  are  repeated  for  different  values  of  second  independent  variable  until  all 
second  independent  variales  have  been  exhausted.  The  next  card  has  an  identi¬ 
fier  corresponds  to  the  second  independent  variable  and  new  values  of  that 
variable  for  the  plane  of  the  second  value  of  third  independent  variable.  The 
values  of  the  first  independent  variable  need  not  be  repeated  if  they  are  the 
same  along  each  line  of  constant  second  independent  variable.  In  each  instance 
where  the  values  are  changed  a  new  card  is  required. 

Last  Table  input  termination  Indicator,  EOT  AS 

All  remaining  tables  for  this  input  section  follow  the  same  pattern  as  above. 

To  end  the  table  read-in  mode,  a  blank  table  reference  number  is  input  be¬ 
hind  the  last  table  of  the  entire  table  set. 


A-3 


NADC-83030-60 


TABLE  A-II.  ONE  PARAMETER  CARD  INPUTS 

COLUMN  LOCATION 

123456789012345678901234567890123456789012345678901234567890 


101 

DRAG 

COEFFICIENT 

VS  MACH 

NUMBER 

z 

1 

0.0 

Y 

1 

0.0 

MACH 

4 

0.0 

0.1 

0.2 

0.3 

CD 

4 

0.010 

0.011 

0.0112 

0.0115 

EOT 
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TABLE  A-III.  TWO  PARAMETER  CARD  INPUTS 

COLUMN  LOCATION 

123456789012345678901234567890123456789012345678901234567890 


104 

Z 

1 

DRAG 

0.0 

COEFFICIENT 

VS  M  AND 

CL 

CL 

3 

0.0 

0.2 

0.3 

MACH 

4 

0.0 

0.1 

0.2 

0.3 

CD 

4 

0.01 

0.02 

0.03 

0.04 

MACH 

3 

0.0 

0.15 

0.2 

CD 

3 

0.1 

0.02 

0.03 

MACH 

3 

0.0 

0.15 

0.2 

CD 

EOT 

3 

0.02 

0.03 

0.04 
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TABLE  A-IV.  THREE  PARAMETER  CARD  INPUTS 


COLUMN  LOCATION 

123456789012345678901234567890123456789012345678901234567890 


226 

CD  VS  M, 

CL,  AND 

CG 

CG 

3 

0.0 

10.0 

20.0 

CL 

3 

0.0 

0.1 

0.2 

MACH 

4 

0.0 

0.1 

0.2 

0.3 

CD 

4 

0.01 

0.01 

0.02 

0.022 

CD 

4 

0.02 

0.02 

0.03 

0.035 

MACH 

3 

0.0 

0.1 

0.3 

CD 

3 

0.03 

0.031 

0.033 

CL 

4 

0.0 

0.15 

0.20 

0.3 

MACH 

3 

0.0 

0.2 

0.3 

CD 

3 

0.011 

0.011 

0.021 

CD 

3 

0.015 

0.015 

0.026 

CD 

3 

0.020 

0.020 

0.036 

CD 

3 

0.025 

0.025 

0.041 

CL 

3 

0.2 

0.3 

0.4 

MACH 

4 

0.0 

0.2 

0.4 

0.6 

CD 

4 

0.01 

0.01 

0.015 

0.020 

MACH 

3 

0.0 

0.2 

0.3 

CD 

3 

0.011 

0.011 

0.022 

CD 

3 

0.021 

0.022 

0.032 

EOT 
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logging  into  the  host  system  and  the  baud  rate  has  been  selected,  different  prompts  will  appear 
depending  on  the  user  response  to  the  initial  interactive  query.  The  response  will  depend  on 
whether  the  user  intends  to  correct  an  existing  file  or  create  a  new  file  via  the  screen  or  tablet. 
After  the  baud  rate  selection  the  next  query  to  appear  will  be: 

(a)  “IS  THIS  A  CREATION  RUN?" 

A  "Y"  response  indicates  a  T APE  1  file  is  to  be  newly  created  and  the  following  prompts  will 
appear. 

(b)  “ENTER  TABLE  TITLE  CARD 

(COLUMNS  1-5  SHOULD  BE  TABLE  REFERENCE  NUMBER)" 

The  use  should  refer  to  the  instructions  in  section  A-2  Table  A-l,  card  1. 

(c)  “ENTER  4  CHARACTERS  FOR  EACH  LABEL  FOR  Z,  Y,  X,  FXYZ 
(separated  by  commas)" 

The  user  should  refer  to  instructions  in  section  A.2  Table  A-l,  card  2. 

(d)  “ENTER  NUMBER  OF  Z  VALUES" 

The  user  should  refer  to  instructions  in  section  A.2  Table  A-l,  card  2. 

(e)  “ENTER  Z  VALUES  IN  ASCENDING  ORDER" 

The  user  should  now  enter  the  values  of  the  Z  parameter  with  blanks  or  commas  between  the  data 
pieces. 

(f)  “WANT  TO  SPECIFY  DECIMAL  PLACES  ON  TAPE2?" 

An  “N"  response  will  by  default,  set  the  number  of  places  at  the  maximum  allowable.  If  the  user 
enters  a  "Y",  this  will  be  prompt  query  (g). 

(g)  “ENTER  NUMBER  OF  DECIMAL  PLACES  FOR  Z,Y,X,FXYZ" 

The  user  should  specify  the  number  of  decimal  places  (up  to  9)  separated  by  blanks  or  commas  for 
Z,Y,X,FXYZ  parameters. 

At  this  point,  the  following  message  will  appear: 

“NO  DATA  TO  BE  FOUND.  .  .ENTER  COMMAND" 

The  user  may  now  begin  creating  the  tabular  data  set  with  either  a  “N"  (new  line)  command 
or  a  “T"  (tablet  operation)  command.  The  reader  is  referred  to  sections  A.5  and  A.6  for  additional 
information. 

If  the  response  to  the  initial  query  (query  (a)  above)  is  “N"  then  this  means  that  data  on  the 
TAPE1  is  to  be  used  and  queries  (b)  thru  (e)  are  skipped. 

A.6  TABLET  INITIATION  PROCEDURE 

Tablet  commands  are  issued  using  a  command  menu.  The  command  menu  is  a  section  of  the 
tablet,  20  one  inch  squares  (10  columns  by  2  rows)  in  which  keying  the  coordinates  within  a 
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square  is  interpreted  as  the  indicated  command.  Figure  A-4  shows  the  positions  of  the  menu 
commands  within  the  10  by  2  inch  squares. 

Tablet  operation  begins  initially  by  attaching  the  command  menu  at  any  convenient  location 
on  the  tablet.  The  menu  should  be  approximately  parallel  with  the  lower  edge  of  the  tablet.  Upon 
first  entering  the  ‘‘T"  command,  the  user  enters  the  position  of  the  upper  left  corner  of  the  menu. 

The  coordinates  of  this  position  are  then  used  to  determine  the  commands.  Next,  the  user  attaches 
the  graph  that  is  to  be  used  at  a  convenient  tablet  location.  Squaring  the  graph  is  not  necessary 
since  any  angular  correction  required  is  performed  in  TIGS.  The  user  then  may  select  any  conven¬ 
ient  orthogonal  axes  and  enters  the  following  requested  information: 

a)  the  coordinate  position  of  the  crossing  point  of  the  orthogonal  axes,  and  the  coordinate 
values  X  and  FXYZ  respectively  at  the  crossing  point. 

b)  the  coordinate  position  of  any  X  axis  point  and  its  value,  (usually  this  coordinate  posi¬ 
tion  is  the  maximum  axis  length) 

After  these  entries  have  been  made  the  "NO  DATA  FOUND  TO  PLOT"  message  will  appear. 
At  this  time  the  user  may  issue  commands  from  the  tablet  menu. 

A.7  TIGS  INTERACTIVE  COMMANDS 

Commands  from  the  screen  involve  only  positioning  the  cross  hairs  and  keying  the  appropriate 
command.  Commands  from  the  tablet  involve  first  selecting  the  command  from  the  menu  and 
then  indicating  the  coordinate  position.  Once  a  tablet  command  has  been  set  it  remains  settjntii 
changed.  The  user  is  free  to  change  to  and  from  tablet  and  screen  command  modes.  The  following 
commands  are  available: 

"A"-  add  point  after.  The  user  positions  the  cross  hairs  and  keys  the  "A"  command  (or 
indicates  the  tablet  command  and  position).  The  system  will  respond  by  drawing 
the  symbol  at  the  new  point.  (Note.  See  the  "C"  command  for  further  discussion.) 

"B"-  add  point  before.  This  command  is  exactly  like  the  "A"  command  except  that  the 
point  is  added  before  the  pointer  position. 

"C"-  position  the  pointer  to  the  array  location  that  the  user  wishes  to  add  a  new  point. 

The  next  command  following  the  "C"  command  to  add  a  point  may  be  an  "A"  to 
add  after  or  a  "B"  to  add  before  the  pointer  position.  In  addition  an  "M"  command 
may  be  used  to  move  to  a  new  location  the  point  indicated  by  the  position  pointer. 

A  "V"  command  may  also  be  used.  It  should  be  noted  that  the  pointer  position  after 
each  added  point  becomes  the  position  of  the  added  point.  Possible  valid  commands 
would  be  "CAABBAVVVAA"  permitting  the  user  to  continuously  add  new  points 
very  rapidly.  Any  other  command  drops  the  pointer  position,  which  must  be  restored 
by  another  "C"  command  to  add  new  points. 

"D"-  delete  the  point  closest  to  the  cross  hairs  or  pen  position. 

"E"-  end  or  terminate  this  plot  and  return  to  TIGS  for  next  plot  if  any. 

"F"-  format  or  change  type  of  curve  drawn  for  each  line  as  follows: 
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ABCDEGHIMP 
add  add  position  delete  end  grid  halt  initial  new  plot 
after  before  pointer  point  plot  switch  tablet  tablet  line  data 

R  S  V  W 

restore  show  value  window 
window  value  input  data 


FIGURE  A-4.  TABLET  COMMAND  MENU 
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ITiP-  switch  determining  type  of  curve  drawn  (iTIP=2,  default)  ITIP 
can  have  the  following  values:  (Note:  a  negative  value  will  have  the 
same  meaning  except  no  symbols  are  drawn.) 

0  indicates  symbols  only,  no  curve  drawn 

1  indicates  linear  fit 

2  indicates  smooth  spline  like  fit  with  respect  to  x  axis 

3  same  as  ITIP=2  except  with  respect  to  y  axis 

4  indicates  data  is  multivalued  and  the  fit  is  with  respect  to  arc  length 
along  curve 

5  indicates  data  is  multivalued  and  forms  a  closed  figure;  the  fit  is  with 
respect  to  arc  length  and  joined  at  the  ends. 

NOTE:  if  data  is  not  in  ascending  order  when  ITIP  =  2  or  ITIP  =  3,  then  curve  fit 
will  default  to  ITIP  =  4;  DEFAULT  format  is  ITIP  =  2. 

"G”-  IGRID  switch-  turns  grid  from  on  to  off  or  off  to  on. 

"H”-  halt  tablet  and  returns  control  to  screen.  This  only  applies  to  tablet  modes. 

"I"-  initialize  tablet  starting  with  graph  coordinate  locations.  This  only  applies  to 
tablet  mode. 

"M"-  move  the  point  indicated  by  the  pointer  position  to  the  new  coordinates  indicated 
by  the  cross  hairs. 

"N"-  begin  a  new  curve  at  point  indicated.  Whenever  this  command  is  issued,  the  user 
will  also  enter  the  new  curve  value.  The  pointer  position  becomes  the  new  point 
permitting  commands  such  as  "NAAMAAAABBBVAB". 

"P"-  re-plot  data  with  scale  as  shown,  see  Note 

"R"-  rescale  data  to  the  largest  size  and  re-plot,  see  Note 

"S"-  show  the  current  coordinate  values  at  the  position  indicated. 

"V"-  values  input;  same  as  "A"  or  "B"  command  except  the  actual  X  and  F  coordinate 
values  are  input. 

"W”-  window  the  plot  to  fit  within  the  rectangle  indicated  by  the  diagonal  between  two 
coordinate  positions  (two  points  are  sent). 

Note:  The  commands  "F",  "P"  and  "R"  have  a  dual  meaning.  If  the  vertical  cross  hair 
is  to  the  left  of  the  vertical  plot  axis  then: 

"F"  indicates  to  change  the  ITIP  of  the  curve  indicated  by  the  position  of  the  horizontal 
cross  hairs.  "P"  indicates  activate  TIGS  to  plot  the  curve  indicated  by  the  position  of  the 
horizontal  cross  hairs.  This  action  is  cumulative  in  that  one  or  more  of  many  curves  may  be 
indicated  in  this  manner.  "R"  is  used  to  deactivate  this  special  mode  and  thus  restore  all  curves  to 
plotting  status. 
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A.8  EXAMPLE  INTERACTIVE  SESSION 

In  this  example  the  TIGS  system  is  used  to  create  a  plot.  The  figures  in  this  example  are 
actual  copies  of  what  the  user  would  see  on  the  Tektronix  screen.  In  these  figures  a  "?"  followed 
by  data  indicates  that  these  data  were  the  user's  response  to  the  indicated  query.  In  the  following 
discussion  numbers  enclosed  by  circles  refer  to  corresponding  numbers  on  a  figure  pointing  to  a 
feature  under  discussion.  In  figure  A-5  (T)  is  the  command  used  to  begin  execution  of  TIGS. 

(2)  is  the  user  response  to  the  query  as  to  transmission  line  rate.  @  indicates  that  a  creation 
is  requested  and  results  in  queries  (4)  thru  @  .  (4)  is  the  main  title  on  the  plot  preceded  by 
the  table  reference  number.  The  table  reference  number  should  be  a  5  digit  integer  number.  The 
title  can  be  up  to  4  lines.  The  user  may  indicate  a  new  line  by  leaving  three  consecutive  blanks 
between  words.  (5)  is  the  response  to  the  query  requesting  four  variable  names  for  the  respective 
data.  Note  that  each  variable  name  must  be  4  characters  in  length;  blanks  count  as  characters.  @ 
is  the  response  to  the  number  of  Z  variables  requested.  Each  Z  value  represents  a  single  plane. 
is  the  response  to  input  each  Z  value.  @  is  the  response  to  the  decimal  place  query  related  to  the 
TAPE?  file.  This  file  is  an  output  file  containing  all  of  the  data  generated  during  this  session.  Each 
prompt,  as  shown,  indicates  that  the  data  on  TAPE?  will  contain  the  maximum  decimal  places  that 
will  fit  with  each  space.  Some  caution  is  necessary  if  the  user  specifies  the  number  of  decimal 
places  for  each  parameter;  precision  could  be  lost  if  a  low  number  of  decimal  places  is  initially 
selected.  A  good  technique  is  to  examine  the  TAPE?  file  with  the  maximum  decimal  places  speci¬ 
fied  first  and  then  re-enter  TIGS,  if  necessary,  and  specify  decimal  places  as  required. 

When  the  user  responds  to  (s)  ,  figure  A-6  will  be  displayed.  The  meaning  of  figure  A-6  is 
that  a  plot  command  as  implied  and  that  data  was  not  found  to  plot.  This  is  a  proper  response  since 
the  user,  through  (5)  ,  on  figure  A-5,  elected  a  creation  run  and  there  is  no  data  as  yet  to  plot. 

The  user  will  notice  for  the  first  time  that  cross  hairs  also  have  appeared  on  the  screen.  (Note  these 
are  not  shown  in  figure  A-6).  The  cross  hair  is  a  prompt  signal  that  an  input  is  reque«ted.  The  input 
is  a  single  upper  case  letter.  A  "RETURN"  is  not  required  after  typing  the  single  letter  command. 
The  single  letter  command  issued  in  this  example  was  a  "N"  indicating  a  new  line.  The  response 
to  this  command  shown  in  figure  A-?.  The  first  prompt,  @  ,  requests  one  set  of  data  coor^ates, 
X  and  FXYZ,  for  one  point.  The  purpose  of  this  is  to  scale  the  final  plot.  The  response  at  is 
a  value  assigned  to  this  line  (this  one  set  of  coordinates  is  the  beginning  of  a  potential  curve). 

After  entering  the  number  one  for  this  query  as  noted  by  ,  the  screen  will  appear  as 

shown  in  figure  A-8.  The  "Y"  shown  on  the  left  top  of  this  figure  is  the  4  character  label  entered  in 

figure  A-5.  The  "^nd  number  under  the  "Y"  is  the  symbol  for  the  first  line  and  the  line 
value  assigned  to  it  as,  in  figure  ?.  The  "A"  in  the  center  of  the  plot  at  (0,0)  is  the  first  point 
(and  only  point)  of  line  aT  The  pointer  positioned  message  indicates  that  the  reference  point  from 
which  to  add  points  has  been  identified.  This  occurred  automatically  since  only  one  point  at  this 
time  is  in  the  plot,  the  first  point.  All  of  the  other  4  character  labels,  including  the  main  plot  label, 
are  also  shown.  The  value  of  the  plot  plane  (Z  value)  is  zero  and  is  shown  at  the  top  right. 

In  the  next  steps  the  user  has  moved  the  cross  hairs  and  "keyed"  the  "A"  characters  indicating 
"add  point  after".  The  "add  point  after"  in  this  context  means  that  the  data  point  storage  of  the 
new  point  is  after  the  point  indicated  by  pointer  position.  The  curves  are  always  drawn  in  the  order 
towards  the  "after"  point.  After  each  point  is  added  the  pointer  position  becomes  the  position  of 
the  added  point.  Figure  A-9  shows  the  addition  of  added  four  points  as  they  would  appear  on  the 
screen.  Figure  A-10  is  a  replot  of  the  data  resulting  from  the  user  keying  a  "P".  This  command 
simply  plots  a  curve  through  the  data  points  shown. 

Figure  A-1 1  is  a  resize  and  replot  resulting  from  an  "R"  command.  The  plot  axes  have  been 
rescaled  to  permit  the  largest  plot  of  the  data  points  that  will  fit  within  the  screen. 
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TIGS 

ENTER  BAUD  RATE  CODE 
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ENTER  TABLE  TITLE  CARO 
(COLUMNS  1-5  SHOULD  BE  TH£  TABLE  REFERENCE  NUMBER) 
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ENTER  4  CHARACTERS  FOR  EACH  LABEL  FOR  Z,Y.X,FXYZ 
(SEPARATED  BY  COMMAS) 
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ENTER  NUMBER  OF  Z  VARUBLES - FREE  FORM 
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TO  SPECIFY  DECIMAL  PLACES  ON  TAPE77 
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FIGURE  A-5.  EXAMPLE  PLOT 
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NO  DATA  POUND  TO  PLOT  ..ENTER  CONNMD 


FIGURE  A-6.  EXAMPLE  PLOT 
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FIGURE  A-7.  EXAMPLE  PLOT 
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FIGURE  A-10.  EXAMPLE  PLOT 
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In  figure  A-1 2  the  cursor  was  first  positioned  at  approximate  co-ordinate  locations  of  .1 5  and 
.16  and  a  "C"  command  issued.  The  "C"  indicates  to  identify  the  closest  point  to  the  intersection 
of  the  cross  hairs  as  pointer  position.  This  command  also  resulted  in  the  message  stating  "pointer 
positioned"  at  the  top  left.  The  four  points  shown  on  figure  12  were  then  added  by  the  user  moving 
the  cursor  and  "Keying"  the  "A"  command.  In  figure  A-13  a  "P"  command  was  issued  first  and 
then  the  cross  hairs  were  at  the  position  indicated  by  the  "B"  symbol  and  the  "N"  command  was 
keyed.  This  resulted  in  the  "input  line  value"  query  shown  in  the  top  left  of  figure  A-13.  In  figure 
A-1 4  the  query  response  is  shown  and  the  user  has  inputted  more  points  by  just  moving  the  cursor 
and  using  the  "A"  command. 

In  figure  A-1 5  the  user  has  replotted  the  data  with  a  "P"  command  and  then  the  cursor  was 
positioned  near  the  end  "B"  point  at  X=.14  and  a  "C"  command  was  keyed.  The  next  command 
sent  by  the  user  was  a  "V",  to  input  an  exact  value.  This  prompted  the  query  "Input  X,Y"  to 
appear.  The  last  query  "A  or  B  mode?"  simply  request  that  the  user  identify  where  in  the  data 
storage  is  the  new  data  point  stored,  before  the  pointer  or  after  the  pointer. 

Figure  A-16  is  a  final  plot  of  the  data  showing  the  new  point.  At  this  point  the  user  keyed  an 
"E"  command  and  "ended"  the  execution.  In  figure  A-1 7,  the  output  file  created  during  this 
example  is  listed  using  the  CED  text  editor,  showing  all  the  data  points. 
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FIGURE  A-12.  EXAMPLE  PLOT 
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FIGURE  A-14.  EXAMPLE  PLOT 
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FIGURE  A-15.  EXAMPLE  PLOT 
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FIGURE  A-17.  EXAMPLE  PLOT 
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CTIG3 

PROGRAM  TIGS (INPUTr 101 , OUTPUT, TAPE1=101 ,TAPE7=101 ,TAPE5=INPUT) 
^  ***»»»***»»«»**»»*»«*****»»*»»»»*«»* 

tigs  tflot  interactice  graphics  system 

Q**** 

C****  M  CADDY  JAN  30  78 

DIMENSION  LT(7) ,XV(30) ,NPTS(30) ,X(300) ,Y(300) ,Z(30) ,A(99) 

DATA  NPLOT/0/ 

DATA  NT/8/ 

DATA  ITIP,IGRID/2,1  / 

DATA  NPTS,XV/30*0,30*0./ 

DATA  IEND/10HE0T  / 

10  FORMAT(A5,7A10) 

REWIND  1 
REWIND  7 
PRINT  20 

20  FORMATC*  TIGS  VER  2.0  8/2/78  *, 

1  /*  IF  THIS  IS  A  CREATION  RUN  ENTER  Y*) 

READ  30, IC 

30  FORMAT(IRI) 

IC=IC-30B 

IF(IC.EQ.I)  GO  TO  40 

C»***  FILE  IS  NOT  BEING  CREATED  READ  IT  FROM  TAPE1 

Q**** 

31  READ  (1,10)  LNO,LT 

IFCLNO.EQ.lOH  )  GO  TO  251 

Ca***»  CALL  IN  Z  VALUES 

CALL  TABR(LZ,NZ,Z,1) 

(;»*»» 

C**»*  READ  IN  Y,X,FXYZ  DATA  FOR  NON  CREATION  RUN 

CALL  TABR(LY,NY,A,1) 

CALL  TABR(LX,N,X,1) 

CALL  TABR(LF,N,Y,1) 

GO  TO  100 
40  PRINT  50 

50  FORMAT(»  ENTER  TABLE  TITLE  CARD*/, 

1  *  (COLUMNS  1-5  SHOULD  BE  THE  TABLE  REFERENCE  NUMBER)*) 

READ  10,LN0,LT 

Q***» 

C****  READ  TABLE  NUMBER  AND  TITLE 

PRINT  60 

60  FORMAT (*  ENTER  4  CHARACTERS  FOR  EACH  LABEL  FOR  Z,Y,X,FXYY*/ 


TIGS0002 

TIGS0003 

TIGS0004 

TIGS0005 

TIGS0006 

TIGS0007 

TIGS0008 

TIGS0009 

TIGS0010 

TIGS0011 

TIGS0012 

TIGS0013 

TIGS0014 

TIGS0015 

TIGS0016 

TIGS0017 

TIGS0018 

TIGS0019 

TIGS0020 

TIGS0021 

TIGS0022 

TIGS0023 

TIGS0024 

TIGS0025 

TIGS0026 

TIGS0027 

TIGS0028 

TIGS0029 

TIGS0030 

TIGS0031 

TIGS0032 

TIGS0033 

TIGS0034 

TIGS0035 

TIGS0036 

TIGS0037 

TIGS0038 

TIGS0039 

TIGS0040 

TIGS0041 

TIGS0042 

TIGS0043 

TIGS0044 

TIGS0045 

TIGS0046 

TIGS0047 
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1  »  (SEPARATED  BY  COMMAS)*) 

C**** 

C***»  READ  TITLES  FOR  EACH  VARIABLE  4  CHARACTERS  LONG 

READ  70,LZ,LY,LX,LF 
70  F0RMAT(4(A4,1X)) 

Q**** 

C****  GET  NUMBER  OF  Z  VARIABLES  AND  VALUES 

C»*** 

PRINT  80, LZ 

80  FORMATC*  ENTER  NUMBER  OF  *,A4,»  VARIABLES -  FREE  FORM*) 

CALL  GETIN(1,Z) 

NZrZd) 

PRINT  90, LZ 

90  FORMATC*  ENTER  *,A4,*  VALUES  , ASCENDING  ORDER— FREE  FORM*) 
CALL  GETIN(NZ,Z) 

C****  WRITE  TO  TAPE7  TITLE  CARD  AND  TABLE  NUMBER 

100  WRITE(7,10)  LNO,LT 
PRINT  110,LZ,LY,LX,LF 

110  FORMATC*  ENTER  NUMBER  OF  DECIMAL  PLACES  FOR  *,4(A4,1X) 

1  ,*  FREE  FORM*) 

C****  GET  NUMBER  OF  DECIMAL  PLACES  FOR  EACH  VARIABLE 

Q**** 

CALL  GETIN(4,XV) 

LZDP=XVC1) 

LYDPrXV(2) 

LXDP=XV(3) 

LFDP=XV(4) 

C****  WRITE  TO  TAPE7  THE  Z  VALUES  ETC... 

Q»*** 

CALL  TF0RM(1,LZ,NZ,Z,LZDP,7) 

c»***  INIALIZE  TEK  SOFTWARE 

CALL  INITT(120) 

CALL  TERMC3,4096) 

CALL  CHRSIZ(4) 

DO  250  IZ=1,NZ 

c****  IP  CREATION  MODE  THEN  SET  DEFAULTS  TO  0 

C***« 

IFCIC.NE.I)  GO  TO  120 

NPTSCDrO 

X(1)=0. 
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Y(1)=0. 

GO  TO  210 

Q**** 

c****  non  creation  mode 

120  CONTINUE 
K=1 

C****  TRANSFER  SECOND  INDEPENDENT  VARIABLE  TO  XV  ARRAY 

Q**** 

DO  130  J=1,NY 
130  XV(J)=A(J) 

IF(IZ.EQ.I)  GO  TO  140 
CALL  TABR(LX,N,X,1) 

CALL  TABR(LF,N,Y,1) 

140  LNX=N 
LNY=N 
NPTS(1)=N 
NPTS(2)=0 

C****  READ  NEXT  SET 

Q**** 

150  CALL  TABR(LW,N,A,1) 
c»»»*  CHECK  FOR  NEXT  Z  GROUP 

C***» 

IF(LW.EQ.LY)  GO  TO  210 

C»»»* 

C**»*  check  for  END  OF  TABLE 

C»»*» 

IF(LW.EQ.4HE0T  )  GO  TO  210 

C**** 

C***»  CHECK  FOR  NEXT  X  DATA 

C»*** 

IF(LW.NE.LX)  GO  TO  170 

C*»*» 

c*»*»  data  is  X  DATA  STORE  IT 

C*»*» 

LOX=LNX 
DO  160  J=1,N 
LNX=LNX+1 
160  X(LNX)=A(J) 

GO  TO  150 

C»»»» 

C**»»  data  had  better  BE  LY 

c»«»* 

170  IF(LW.NE.LF)  STOP 
c**»» 
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C****  IF  data  has  not  be  input  for  X  DATA  USE  LAST  VALUES 

Q**** 

IF(LNX.GT.LNY)  GO  TO  190 

LLrLOX 

DO  180  Jr1,N 

LNXrLNX+1 

LL=LL+1 

180  X(LNX)=X(LL) 

Q**** 

C»**»  UPDATE  COUNTERS 

Q**** 

190  K=K+1 
NPTS(K)=N 
NPTS(K+1)=0 

Q**** 

C****  LOAD  Y  DATA 

DO  200  J=1,N 
LNY=LNY+1 
200  Y(LNY)=A(J) 

C****  GO  BACK  TO  GET  NEXT  GROUP 

C**** 

GO  TO  150 

c***»  plot  data 

c«*** 

210  CALL  TIGPPR(NPLOT,LF, 1,LX, 1,LT,8,X,Y,NPTS,LY, 1,XV,LYDP,ITIP, 
1  IGRID,LZ,Z(IZ)) 

CALL  ANMODE 

c**»*  COUNT  NUMBER  OF  Y  VALUES 

C***» 

NY=0 

DO  220  1=1,30 

IF(NPTS(I).EQ.O)  GO  TO  230 
NY=NY+1 
220  CONTINUE 
GO  TO  250 

C»»»* 

C»**»  WRITE  TO  TAPE?  Y  DATA  ETC.... 

C»*** 

230  CALL  TFORM(1,LY,NY,XV,LYDP,7) 

LOC=1 

J=0 

2^10  J=J+1 

NP=NPTS(J) 

IF(NP.EQ.O)  GO  TO  250 
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Q**** 

c*«**  miE  TO  TAPE?  X  DATA  ETC... 

Q**** 

CALL  TF0RM(L0C,LX,NP,X,LXDP,7) 

Q**** 

C»*»*  WRITE  TO  TAPE?  Y  DATA  ETC... 

CALL  TFORM(LOC,LF,NP,Y,LFDP,?) 

LOC=LOC+NP 
GO  TO  240 

250  CONTINUE 

WRITE  (?,10)  lEND 

C»***  IF  MON  CREATION  MODE  THEN  GO  BACK  TO  READ  NEXT  TITLE 

IF(IC.NE.I)  GO  TO  31 

251  WRITE  (?,10) 

RE’^IND  ? 

END 

CTFORM 

SUBROUTINE  TFORM(LOC,UB,N,X,  IP,K) 

DIMENSION  X(1),IFORM(3) 

FORMATTING  SUBROUTINE  FOR  TPLOT  FORMAT 
LOC  IS  THE  LOCAL  ARRAY  POSITION  TO  PRINT  FRCM 
C*»*»  lab  is  the  4  CHARACTER  LABEL 

C****  N  IS  THE  NUMBER  TO  PRINT 

C**»»  X  IS  THE  ARRAY  CONTAINING  THE  VALUES 

IP  IS  THE  NUMBER  OF  DECIMAL  PUCES  TO  USE  IN  FORMAT 
IF(IP.LT.O)  IP=0 
IFCIP.GT.g)  IP=9 
JO=LOC-1 
NP=N 

IFCNP.GT.?)  NP=? 

IF0RM(1)=10H(A4,I3,3X, 

IF0RM(2)=55555555420634335?33B+IP 

IFORM(3)=10H) 

WRITE (K , IFORM )  LAB , N , ( X ( 1+ JO ) , 1= 1 , NP ) 

IFORM(1)=10H(  10X, 

IFCN.GT.?)  WRITE (K, IFORM )(X(I+JO), 1=8, N) 

RETURN 

END 

CTABR 

SUBROUTINE  TABR(LAB,N,A,K) 

DIMENSION  A( 1 ) 

READ(K,10)  LAB,N,(A(I),I=1,?) 

10  FORMAT(A4,I3,3X,?F10.0) 

IFCN.GT.?)  READ(K,20)  (A(I),I=8,N) 

20  F0RMAT(10X,?F10.0) 
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30  IS^JB(I)=I-^1 

MERGE  HERE  TO  REPLOT 

40  CALL  BINITT 
LCNT=3120 

IGRID1r(3»IGRID+7)*.5 

SUM  UP  NUMBER  OF  POINTS 

NL=0 
NPT0T=0 
DO  60  1=1,30 
NrNPTAd) 

IF(N.EQ.O)  GO  TO  70 
NL=NL+1 

60  NPTOTrNPTOT-i-N 

SET  STORAGE  LIMIT  TO  NPTOT  FIRST  PASS 
70  IF(NSTOR.EQ.O)  NSTOR=NPTOT 
IF(NPTOT.GT.O)GO  TO  90 
NSTORrO 

CALL  MOVABS(0,LCNT) 

CALL  ANMODE 
PRINT  80 

80  FORMATC*  NO  DATA  FOUND  TO  PLOT  ..ENTER  COMMAND*) 
LCNT=LCNT-LDEL 
IPL0T=0 
GO  TO  200 


SECOND  INDEPENDENT  VARIABLE  TITLE 

90  IF(NCC.LE.O)GO  TO  140 
CALL  MOVABS(0,LCNT) 

CALL  ANMODE 

PRINT  110,(LABVAL(J1),J1=1,NCC) 
CALL  MOVABSC 2800, 2800) 

CALL  ANMODE 
PRINT  100,LZ,ZVAL 
100  FORMAT(A4,*=*,G13.5) 

110  FORMAT(8A10) 

LCNT=LCNT-LDEL 

KL=0 

KH=55B 

DO  130  J1=1,NL 
LCNT=LCNT-LDEL 
CALL  MOVABSC 0,LCNT) 

KL=KL-t-1 
CALL  ANMODE 

PRINT  120,KH,KL,VLABL(J1) 


TIGP0047 

TIGP0048 
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TIGP0075 

TIGP0076 

TIGP0077 

TIGP0078 

TIGP0079 

TIGP0080 

TIGP0081 

TIGP0082 

TIGP0083 

TIGP0084 

TIGP0085 

TIGP0086 

TIGP0087 

TIGP0088 

TIGP0089 

TIGP0090 

TIGP0091 

TIGP0092 

TIGP0093 

TIGP0094 

TIGP0095 

TIGP0096 

TIGP0097 

TIGP0098 

TIGP0099 
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RETURN 

END 


CTIGP 


TEK  INTERACTIVE  GPPR  M  CADDY  FEB  78 


SUBROUTINE  TIGPPR(NPL0T,LABY,N1,LABX,N2,LABTL,NT  ,X,Y, 

1  NPTA , LABVAL , NCC , VLABL , NDECVIN , ITIP , IGRID , LZ , ZVAL) 
COMMON/TKTRNX/ITEKC(60 ) 

DIMENSION  X(200 ) , Y( 200 ) , LABTL  C 9 ) , NPTA( 30 ) , VLABL ( 30 ) , 

1  LABX(5) ,LABY(5) ,VTEM(8) ,LABVAL(8) , IQUICK(30) ,ISUB(300) 

DIMENS ION  MSG 1 ( 20 ) , MSG2 ( 20 ) , MSG4 ( 1 0 ) , MSG5 ( 1 0 ) , MSG6 (15), lALTM (6,2) 
EQUIVALENCE  ( BEG ( 1 ) , XBEG ) , ( BEG ( 2 ) , YBEG ) 

EQUIVALENCE  (DEL(1 ),DELX) ,(DEL(2) ,DELY) ,(ITAB,LTV(2)) 

EQUIVALENCE  ( EN ( 1) , XEND ) , ( EN ( 2 ) , YEND ) 

EQU I VALENCE  ( lOFF , ITEKC ( 30 ) ) , ( TXMIN , ITEKC ( 1 ) ) 

C  SET  LINE  SPACING 

COMMON/TEKGPPR/LDEL , LCNT , MAXSR , LTV ( 1 7 ) , EN ( 2 ) , DEL ( 2 ) , BEG ( 2 ) , RDX2 , 

1 RDY2 , NLINE , NDRAW( 30 ) ,MODE ( 30 ) 

DATA  MSG1/46,q46,80,111,105,110,1l6,101,114,32, 

1  80,111,115,105,116,105,111,110,101,100/ 

DATA  MSG2/73, 110, 112, 117, 116, 32, 76, 105, 110, 101, 

1  32,86,97,108,117,101,32,32,32,32/ 

C»»»»  illegal  MESSAGE 

DATA  MSG4/42, 73, 108, 108, 101, 103, 97, 108, 42, 32/ 

DATA  MSG5/73, 110, 112, 117, 116, 32, 88, 44, 89, 32/ 

DATA  MSG6/65, 32, 111, 114, 32, 66, 32, 109, 111, 100, 

1  101,63,32,32,32/ 

DATA  ((IALTM(I,J),I=1,6),J=1,2)/65,102,116,101,114,32, 

1  66,101,102,111,114,101/ 

DATA  IQUICK/0,0,1,1,2,3,10,4,0,0, 

1  0 , 0, 0 , 5 , 0, 6 , 0 , 7, 8 , 0 , 

2  0,0, 9, 0,0, 0,0, 0,0,0/ 

A=1./KIN(1.) 

LDEL=50 

IF(NPLOT.GT.O)  GO  TO  20 
DO  10  1=1,8 
10  LTV(I)=0 

20  IWIN=0 
NPL0T=NPL0T+1 
IGRID=1 

DO  21  1=1,30 

21  MODE(I)=ITIP 
NTL=NT 
NLINE=0 
NST0R=0 

C  SET  STORAGE  POINTER  TO  INITIAL  SEQUENCE 

DO  30  1=1,299 
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120  FORMAT(1X,2R1,G13.5) 

130  CONTINUE 
C 

C  PREPARE  TEKTRONIX  AGII  COMMON 

C 

140  CONTINUE 

IPL0T=IPL0T+1 
CALL  CHRSIZ(4) 

C  SET  SCREEN  WINDOW  SIZE 

CALL  SLIMX( 640, 4000) 

CALL  SLIMYC 300, 2700) 

C  SET  TICK  SIZES 

C  CALL  XTICS(14) 

C  CALL  YTICSdO) 

IF(IWIN.NE.O)  GO  TO  170 

AXMAX=-1.E99 

AYMAX=-1.E99 

AXMIN=+1.E99 

AYMIN=+1.E99 

C  SET  MIN  AND  MAX  DATA  VALUES 

K=1 

DO  150  I=1,NPT0T 
AXMIN =AM IN  1 ( AXMIN , X(K » 

AYMIN  =AMIN 1 ( AYMIN , Y(K ) ) 
AXMAXrAMAXI ( AXMAX, X(K ) ) 

A YMAX=AMAX 1 ( AYMAX , Y (K ) ) 

KLASTrK 

C  SET  KLAST  TO  END  STORAGE  VALUE 

150  K=ISUB(K) 

IWIN=1 

IF ( AXMIN. NE. AXMAX)  GO  TO  160 

AXMIN=AXMIN-.5 

AXMAX=AXMAX+.5 

160  IF( AYMIN. NE. AYMAX)  GO  TO  170 
AYMIN=AYMIN-.5 
AYr'1AX=AYMAX+.  5 

C  SET  VIRTUAL  WINDOW 

170  CALL  DLIMXCAXMIN, AXMAX) 

CALL  DLIMYC AYMIN, AYMAX) 

CALL  XLEN(28) 

CALL  YLEN(28) 

CALL  XFRM(IGRIDI) 

CALL  YFRM(IGRIDI) 

NBASE=IBASEX(0) 

DO  180  Ir1,2 
CALL  LOPTIM(NBASE) 

CALL  WIDTH (NBASE) 

CALL  SPREAD(NBASE) 
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CALL  TSET(NBASE) 

180  NBASE=IBASEY(0) 

EN ( 1) rCOMGET ( IBASEX(27 ) ) 

EN(2)=COMGET(IBASEY(27)) 

BEGd  )=COMGET(IBASEX(26)) 

BEG(2)=COMGET(IBASEY(26)) 

DELXr (XEND-XBEG)/3360 . 

DELY=(YEND-YBEG)/2400. 

C 

C  FIND  VIRTUAL  SPACE  TO  SCREEN  SPACE  SCALING  PARAMETERS 
C 

RDX2=1./(DELX*DELX)  ' 

RDY2=1./(DELY*DELY) 

CALL  SEIWIN 
CALL  GRID 

CALL  LABELdBASEYCO)) 

CALL  LABEL(IBASEX(0)) 

CALL  DRAWIT(NL,NPTA,X,Y,ISUB) 

C 

C  AXIS  LABELS 

C 

CALL  CHRSIZ(3) 

CALL  TTITE( 2320, 3000, NTL,LABTL, 80,0) 

CALL  TTITE(2320, 100, N2,LABX, 80,0) 

CALL  TTITE (450, 1500, N1,LABY, 80,1) 

C 

C  MERGE  HERE  FOR  INTERACTIVE  FUNCTIONS  (BELL) 

C 

200  IF(LCNT.LT.220)  GO  TO  530 
CALL  CHRSIZ(4) 

IF(NPTOT.EQ.I)  GO  TO  240 
CALL  GETVAL ( ICHAR , XO , YO ) 

210  IF(ICHAR.LE.64.0R.ICHAR.GE.95)G0  TO  220 
ICHAR=ICHAR-64 
ICHECK=IQUICK( ICHAR) 

IF(ICHECK.EQ.O)  GO  TO  220 

GO  TO  (300, 400, 440, 200, 460, 500, 540, 560, 590, 455 ),ICHECK 
220  LCNT=LCNT-LDEL 

CALL  NOTATE(0,LCNT,10,MSG4) 

GO  TO  200 
C 

C  ADD  POINT  AFTER  OR  BEFORE  SPECIFIED  POINT  (A  OR  B) 

C 

C  CHECK  IF  C  COMMAND  AND  FIRST  POINT. 

C 

230  IF(NPT0T.EQ.0)G0  TO  460 
240  LCNT=LCNT-LDEL 

CALL  NOTATE (0,LCNT, 20, MSG 1) 
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250  CALL  GETVAL(ICHAR,X0,Y0) 

C  CHECK  FOR  NEW  LINE  COMMAND 

IF(ICHAR.EQ.86)  GO  TO  580 
C  CHECK  FOR  ADD  AFTER 

260  IF(ICHAR.EQ.65)  GO  TO  270 
C  CHECK  FOR  MOVE 

IF(ICHAR.EQ.77)  GO  TO  270 
C  IF  NOT  A  B  OR  M  GO  TO  NEW  COMMAND 

IF(ICHAR,NE.66)  GO  TO  210 
270  CALL  POINTA(XO,YO) 

CALL  MOVEA(XO,YO) 

IF(IOFF.EQ.O)CALL  ANCHO(IS) 

IF(ICHAR.EQ.77)  GO  TO  290 
NPTOT=NPTOT-i-1 

C  INCREMENT  STORAGE  COUNTER 

NSTOR=NSTOR+1 
NPTA  ( ISA  VE )  =NPTA  ( ISA  VE )  -t- 1 
C  MOVE  POINTER  OF  CLOSEST  POINT  TO  END 

ISUB(NSTOR)=ISUB(JSAVE) 

C  CHANGE  CLOSEST  POINTER  TO  ACCESS  LAST  POINT 

ISUB(JSAVE)=NSTOR 
IF(ICHAR.EQ.65)  GO  TO  280 

C  MOVE  OLD  POINT  TO  LAST  POINT  (  INSERT  BEFORE) 

X(NSTOR)=X(JSAVE) 

Y(NSTOR)=Y(JSAVE) 

GO  TO  290 

C  NEW  POINT  ADD  AFTER 

280  IF(KLAST.EQ.JSAVE)  KLAST=NSTOR 
JSAVErNSTOR 
290  X(JSAVE)=X0 
Y(JSAVE)=Y0 
GO  TO  250 

DELETE  POINT  (D) 

300  DSAVE=1.Ei»0 

IF(NPTOT.EQ.O)  GO  TO  200 
IS=64 
NSUMrl 
K=1 

DO  340  1=1, NL 
NEND=NSUM-hNPTA(I)-1 
DO  330  J=NSUM,NEND 
IF(NLINE.EQ.O)  GO  TO  310 
IF(NDRAW(I).EQ.O)  GO  TO  320 
310  XDX=X(K)-X0 
YDY=Y(K)-Y0 

DIST=XDX*XDX*RDX2-t-YDY*YDY*RDY2 
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TIGP0207 

TIGP0208 

TIGP0209 

TIGP0210 

TIGP0211 

TIGP0212 

TIGP0213 

TIGP0214 

TIGP0215 

TIGP0216 

TIGP0217 

TIGP0218 

TIGP0219 

TIGP0220 

TIGP0221 

TIGP0222 

TIGP0223 

TIGP0224 

TIGP0225 

TIGP0226 

TIGP0227 

TIGP0228 

TIGP0229 

TIGP0230 

TIGP0231 

TIGP0232 

TIGP0233 

TIGP0234 

TIGP0235 

TIGP0236 

TIGP0237 

TIGP0238 
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IF(DIST.GE.DSAVE)GO  TO  320 

TIGPO239 

DSAVE=DIST 

TIGP0240 

JSAVErK 

TIGP0241 

ISAVErl 

TIGP0242 

320  KLASTrK 

TIGP0243 

330  K=ISUB(K) 

TIGP0244 

340  NSUM=NEND+1 

TIGP0245 

IS=ISAVE+64 

TIGP0246 

350  IF(IS.LE.90)GO  TO  360 

TIGP0247 

IS=IS-90 

TIGP0248 

GO  TO  350 

TIGP0249 

360  CALL  POINTA(X(JSAVE),Y(JSAVE)) 

TIGP0250 

IF(IOFF.EQ.O)CALL  ANCHO(IS) 

TIGP0251 

IF(ICHAR.NE.4)  GO  TO  230 

TIGP0252 

NPTOT=NPTOT-1 

TIGP0253 

.KzJSAVE 

TIGPO254 

IF  DELETED  POINT  IS  LAST  ONE  SKIP  SHIFT 

TIGP0255 

IF(KLAST.EQ.JSAVE)GO  TO  370 

TIGP0256 

GET  POINTER  OF  NEXT  POINT 

TIGP0257 

K=ISUB(JSAVE) 

TIGP0258 

TRANSFER  POINTER  OF  NEXT  POINT  TO  DELETED  POINT 

TIGP0259 

ISUB(JSAVE)=ISUB(K) 

TIGP0260 

MOVE  VALUE  OF  NEXT  POINT  TO  DELETED  POINT 

TIGP0261 

X(JSAVE)=X(K) 

TIGP0262 

Y(JSAVE)=Y(K) 

TIGP0263 

370  IF(NPTOT.EQ.1)NSTOR=1 

TIGP0264 

ZERO  DELETED  POINTER 

TIGP0265 

ISUB(K)=0 

TIGP0266 

NPTA ( ISA VE ) =NPTA ( ISA VE ) - 1 

TIGP0267 

IF(NPTA(ISAVE).GT.O)GO  TO  390 

TIGP0268 

NPTA ( ISA VE)=0 

TIGP0269 

J=0 

TIGP0270 

DO  380  1=1, NL 

IF(I.EQ.ISAVE)GO  TO  380 

TIGPO27I 

TIGP0272 

J=J+1 

TIGP0273 

NPTA(J)=NPTA(I) 

TIGPO274 

VLABL(J)=VLABL(I) 

TIGP0275 

380  CONTINUE 

TIGP0276 

NPTA(NzL)=0 

TIGP0277 

NL=NL-1 

TIGP0278 

390  GO  TO  200 

TIGP0280 

TIGP0281 

END  (E) 

TIGP0282 

TIGP0283 

400  CALL  NEWPAG 

L=1 

DO  431  I=2,NPTOT 

K=ISUB(L) 

TIGP0284 
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IF(I.EQ.K)  GO  TO  431 
J=K 

JLEFTrNPTOT-i-1-1 
DO  420  KK=1,JLEFT 
IF(J.EQ.I)  GO  TO  430 
JOrJ 

420  JrlSUBCJO) 

430  ISUB(J0)=K 
ISUB(L)=I 
IS=ISUB(I) 
ISUB(I)=ISUB(K) 
ISUB(K)=IS 
XS=X(I) 

X(I)=X(K) 

X(K)=XS 


XS=Y(I)  TIGP0297 

Y(I)=Y(K)  TIGP0298 

Y(K)=XS  TIGP0299 

431  L=I 

RETURN  TIGP0304 

TIGP0305 

FORMAT  (F)  TIGPO3O6 

TIGPO3O7 

440  LCNT=LCNT-LDEL  TIGPO3O8 

CALL  MOVABS(0,LCNT)  TIGP0309 

CALL  ANMODE  TIGPO3IO 

IY=(Y0-YBEG)/DELY-f300  0354 

II=(3045-IY)/50-f1  0355 

IF(II.LT.I)  11=1  0356 

IFCII. GT. ML)  II=NL  0357 

PRINT  450,M0DE(II)  TIGPO3II 

450  FORMATC*  ITIP  =  *,I2)  TIGP0312 

CALL  GETIN(1,VTEM)  TIGPO313 

M0DE(II)=VTEM(1)  TIGP0314 

LCNT=LCNT-LDEL 


IF  F  OUTSIDE  OF  AXIS  THE  SET  ALL  CURVE  MODES 


IF(XO.LE.TXMIN)  GO  TO  200 
DO  451  1=1,30 
451  M0DE(I)=VTEM(1) 

GO  TO  200 


TIGPO317 


CHANGE  GRID  OPTION 


455  IGRID=-IGRID 
GO  TO  200 
C 


TIGPO3I8 
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C  NEW  LINE  (N) 

C 

460  NPT0T=NPT0T+1 
NST0R=NST0R+1 

IF(NPTOT.NE.I)  ISUB(KLAST)=NSTOR 
KLAST=NSTOR 

IF(IPLOT+ITAB.EQ.O)GO  TO  530 

461  X(NSTOR)=XO 
Y(NSTOR)=YO 
NL=NL+1 
IS=NL+64 

C 

C  NO  DATA  THEN  DO[T  SYMBOL  IT 
C 

IF(IPLOT.EQ.O)  GO  TO  490 
470  IF(IS.LE.90)GO  TO  480 
IS=IS-90 
GO  TO  470 

480  CALL  POINTA(XO,YO) 

IF(IOFF.EQ.O)CALL  ANCHO(IS) 

490  NPTA(NL)=1 
NL1=NL+1 
NPTA(NL1)=0 
LCNT=LCNT-LDEL 
CALL  NOTATE (0,LCNT, 20, MSG2) 

LCNT=LCNT-LDEL 
CALL  M0VABS(0,LCNT) 

CALL  ANMODE 

CALL  GETIN(1,VLABL(NL)) 

ISAVErNL 

JSAVE=NSTOR 

IF(IPLOT.EQ.O)  GO  TO  530 
GO  TO  240 
C 

C  PLOT  (P) 

C 

C 

C  CHECK  FOR  TABLET  MODE,  SKIP  SPECIAL  P  SECTION  IF  TABLET 

C 

500  IF(ITAB.EQ.I)  GO  TO  530 
IF(XO.GT.TXMIN)  GO  TO  530 
IF(NLINE.GT.O)  GO  TO  520 
DO  510  1=1, NL 
510  NDRAW(I)=0 
520  IY=(Y0-YBEG)/DELY+300 
II=(3045-IY)/50+1 
IF(II.LT.I)  11=1 
IF(II.GT.NL)  II=NL 


TIGPO319 

TIGPO32O 

TIGPO32I 

TIGPO322 

TIGPO323 

TIGPO324 

TIGPO325 

TIGP0326 

TIGPO327 

TIGP0328 


TIGPO329 

TIGPO33O 

TIGPO33I 

TIGPO332 

TIGPO333 

TIGPO334 

TIGPO335 

TIGPO336 

TIGPO337 

TIGP0338 

TIGPO339 

TIGP0340 

TIGPO341 

TIGP0342 

TIGPO343 

TIGPO344 

TIGPO345 


TIGPO347 

TIGP0348 

TIGPO349 


TIGPO35O 

TIGPO35I 

TIGPO352 

TIGPO353 

TIGPO354 

TIGPO355 

TIGP0356 

TIGPO357 
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NDRAW(II)=1 

TIGP0358 

NLIME=1 

TIGP0359 

GO  TO  200 

TIGP0360 

530  CALL  NE^PAG 

TIGP0361 

GO  TO  40 

TIGP0362 

c 

TIGP0363 

c 

RESTORE  WINDOW  (R) 

TIGP0364 

c 

TIGP0365 

540  IF(XO.GT.TXMIN.OR.NLINE.EQ.O)  GO  TO  550 

TIGP0366 

NLINE=0 

TIGP0367 

GO  TO  200 

TIGP0368 

550  IWINrO 

TIGP0369 

GO  TO  530 

TIGP0370 

c 

TIGP0371 

c 

SHOW  VALUE  (S) 

TIGP0372 

c 

TIGP0373 

560  LCNT=LCNT-LDEL 

TIGP0374 

CALL  MOVABS(0,LCNT) 

TIGP0375 

CALL  ANMODE 

TIGP0376 

PRINT  570,X0,Y0 

TIGP0377 

570  FORMATC  *X=* , G1 3 • 5 , / , *Y=* , G1 3 . 5 ) 

TIGP0378 

LCNT=LCNT-LDEL 

TIGP0379 

GO  TO  200 

TIGP0380 

c 

TIGP0381 

c 

VALUE  IN  (V) 

TIGP0382 

c 

TIGP0383 

580  LCNT=LCNT-LDEL 

TIGP0384 

CALL  NOTATE (0,LCNT, 10, MSG5) 

TIGP0385 

LCNT=LCNT-LDEL 

TIGP0386 

CALL  MOVABS(0,LCNT) 

TIGP0387 

CALL  ANMODE 

TIGP0388 

CALL  GETIN(2,VTEM) 

TIGP0389 

XOrVTEMd  ) 

TIGP0390 

r 

Y0=VTEM(2) 

TIGP0391 

u 

c 

r 

CHECK  FOR  N  COMMAND 

VALUE  INPUT  SECTION. 

Lr 

IF(IPLOT.EQ.O)  GO  TO 

461 

LCNT=LCNT-LDEL 

TIGP0392 

CALL  NOTATE (0,LCNT,1 3, MSG6) 

TIGP0393 

CALL  TINPUT(ICHAR) 

TIGP0394 

GO  TO  260 

TIGP0395 

c 

TIGP0396 

c 

WINDOW  (W) 

TIGP0397 

c 

TIGP0398 

590  CALL  GETVAL(ICHAT,X1 

,Y1) 

TIGP0399 

AXMINrAMIMI  (X0,X1) 

TIGP0400 

AXMAX=AMAX1(X0,X1) 

TIGP0401 

B-15 


N  ADC-83030-60 


AYMIN=AMIN1(Y0,Y1) 

AYMAX=AMAX1(Y0,yi) 

IWIN=1 
GO  TO  530 
END 

CGETVAL 

SUBROUTINE  GETVAL(ICHAR,XV,YV) 

COMMON/TEKGPPR/DUM( 3 ) , ICL , ITAB , HABS, XS , YS , DUM2(20 ) , NLINE 
C  ICL=0  INITIALIZE  TABLET 

C  IC=0  NOT  IN  CONTINUOUS  MODE 

C  ITABrO  SCREEN  CURSER 

C***»  GET  VALUE  AND  CHARACTER  FROM  CROSS  HAIRS  OR  TABLET 

C****  CHECK  FOR  TABLE  INPUTS 

10  IFCITAB.EQ.DGO  TO  20 
CALL  VCURSR(ICHAR,XV,YV) 

C***»  CHECK  TO  SEE  IF  SCREEN  COMMAND  WAS  TO  ACTIVATE  TABLET 

IF(ICHAR.NE.84)GO  TO  30 

ITAB=1 

ICLrITABS 

C****  SET  TABLET  LAST  CHARACTER  (ALSO  USED  AS  A  FLAG  FOR  INITIALIZATION 

20  ICHAR=ICL 

CALL  TABVU(ICHAR,XV,YV) 

C»***  SAVE  LAST  TABLET  CHARACTER  COMMAND 

ICLrICHAR 

NLINE=0 


TIGP0402 

TIGP0403 

TIGP0404 

TIGP0405 

TIGP0406 

GETVOOOl 

GETV0002 

GETV0003 

GETV0004 

GETV0005 

GETV0006 

GETV0007 

GETV0008 

GETV0009 

GETV0010 

GETV0011 

GETV0012 

GETV0013 

GETV0014 

GETV0015 

GETV0016 

GETV0017 

GETV0018 

GETV0019 

GETV0020 

GETV0021 

GETV0022 

GETV0023 

GETV0024 

GETV0025 

GETV0026 

GETV0027 

GETV0028 

GETV0029 

GETV0030 


C 

C  SET  FLAG  TO  PLOT  ALL  LINES  IN  TABLET  MODE 


C 

C»***  GETV0031 

C»***  CHECK  FOR  TABLET  HALT  COMMAND  GETV0032 

C»»**  GETV0033 

IF(ICL.NE.72)GO  TO  30  GETV0034 

C****  GETV0035 

C****  TURN  OFF  TABLET  AND  SAVE  LAST  COMMAND  GETV0036 

C»«*»  GETV0037 

ITAB=0  GETV0038 

ITABSrICL  GETV0039 

RETURN  GETV0040 
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30  IF(ICHAR.EQ.69)ITABS=0 
RETURN 
END 


CTABVU 

SUBROUTINE  TABW(ICHAR,XV, YV) 

COMMON/TEKGPPR/LDEL , LCNT , MAXSR , LTV( 5 ) 

1  ,I^,MX1,MY1,MX2,MY2,XB,YB,FACX,FACY,ANG,MXB,MYB 
DIMENSION  MSG 1(54) ,MSG2(43) ,MSG3(43) ,MSG4(18) ,ICONV(2, 10) , 
1  IRETN(2,10),XTEM(2) 

DATA  ((IC0NV(I,J),J=1, 10), 1=1, 2)765,66,67, 68, 69, 71, 72, 

1  105,78,80,82,83,86,87, 

1  32  32  ,32  32,32  32/ 

DATA  ((IRETN(I,J),J=1,10),I=1,2)/  0,’o,  0,  6,  1,  1,  1, 

1  1,0,  1,  1,  0,  1,  0, 

1  11,1111/ 

DATA  MSG 1/  83,113,117,  97,114,101,  32,109^101,110, 

1  117,  32,119,105,116,104,  32,116,  97,  98, 

1  108,101,116,  32,  97,110,100,  32,116,111, 

1  117,  99,104,  32,117,112,112,101,114,  32, 

1  108,101,102,116,  32,109,101,110,117,  32, 

1  100,111,116,  46/ 

DATA  MSG2/  84,111,117,  99,104,  32,  97,120,105,115, 

1  32,111,114,105,103,105,110,  32,  97,110, 

1  100,  32,101,110,116,101,114,  32,118,  97, 

1  108,117,101,115,  32,  88,  32,  97,110,100, 

1  32,  89,  46/ 

DATA  MSG3/  84,111,117,  99,104,  32,  32,  32,  97,120, 

1  105,115,  32,  97,116,  32,109,  97,120,  32, 

1  108,101,110,103,116,104,  32,  97,110,100, 

1  32,101,110,116,101,114,  32,118,  97,108, 

1  117,101,  46/ 

DATA  MSG4/  76,  97,115,116,  32,  99,111,109,109,  97, 

1  110,100,  32,119,  97,115,  32,  32/ 

IF(ICHAR.NE.O)GO  TO  30 
C  TABLET  HAS  NOT  BEEN  SET  CHECK  IT 

LS=100 

CALL  TABINT(1,0,0) 

CALL  NEWPAG 
LCNT=3120-LDEL 
C  GET  MENU  POSITION 

CALL  NOTATE (0, LCNT, 54, MSG 1) 

CALL  BELL 

CALL  ONEPNT(MX1,MY1) 

MX2=MX1+1000 
MY2=MY1-200 
GO  TO  20 
10  LCNT=3120 
CALL  NE'^PAG 


GETV0041 

GETV0042 

GETV0043 

TABV0001 

TABV0002 

TABV0003 

TABV0004 

TABV0005 

TABV0006 

TABV0007 

TABV0008 

TABV0009 

TABV0010 

TABV0011 

TABV0012 

TABV0013 

TABV0014 

TABV0015 

TABV0016 

TABV0017 

TABV0018 

TABV0019 

TABV0020 

TABV0021 

TABV0022 

TABV0023 

TABV0024 

TABV0025 

TABV0026 

TABV0027 

TABV0028 

TABV0029 

TABV0030 

TABV0031 

TABV0032 

TABV0033 

TABV0034 

TABV0035 

TABV0036 

TABV0037 

TABV0038 

TABV0039 

TABV0040 

TABV0041 

TABV0042 

TABV0043 

TABV0044 

TABV0045 
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20  LCrrfrLCN'T-LDEL 

C  GET  COORDINATE  INTERSECTION 

CALL  NOTATE(0,LCNT,43,MSG2) 

CALL  BELL 

CALL  ONEPNT(MXB,MYB) 

LCNT=LCNT-LDEL 
CALL  MOVABS(0,LCNT) 

CALL  ANMODE 
CALL  GETIN(2,XTEM) 

XB=XTEM( 1 ) 

YB=XTEM(2) 

LCNT=LCNT-LDEL 

MSG3(7)=88 

C  GET  X  AXIS  POSITION  MAX 

CALL  NOTATE(0,LCNT,43,MSG3) 

CALL  BELL 

CALL  ONEPNT(MXI^,NXM) 

LCNT=LCNT-LDEL 
CALL  MOVABS(0,LCNT) 

CALL  ANMODE 

C  GET  VALUE  AT  POSITION 

CALL  GETIN(1,XM) 

DX=MXM-MXB 

DY=NXM-MYB 

C  COMPUTE  ANGLE  CORRECTION 

ANG=ATAN2(DY,DX) 

LCNT=LCNT-LDEL 

MSG3(7)=89 

C  GET  Y  AXIS  POSITION  MAX 

CALL  NOTATE(0,LCNT,43,MSG3) 

CALL  BELL 

CALL  ONEPNT(MYM,NYM) 

LCNT=LCNT-LDEL 
CALL  MOVABS(0,LCNT) 

CALL  ANMODE 

C  GET  VALUE  AT  POSITION 

CALL  GETIN(1,YM) 

DY=NYM-MYB 

COSA=COS(ANG) 

C  SET  UP  COMMON  FACTORS  FOR  ANGLE  CORRECTIONS 

FACX=(XM-XB)»COSA/DX 
FACY=(YM-YB)*COSA/DY 
INITrl 
)0/=XM 
YVrYM 
ICHAR=87 

C  RETURN  PLOT  COMMAND 

RETURN 


TABV0046 

TABV0047 

TABV0048 

TABV0049 

TABV0050 

TABV0051 

TABV0052 

TABV0053 

TABV0054 

TABV0055 

TABV0056 

TABV0057 

TABV0058 

TABV0059 

TABV0060 

TABV0061 

TABV0062 

TABV0063 

TABV0064 

TABV0065 

TABV0066 

TABV0067 

TABV0068 

TABV0069 

TABV0070 

TABV0071 

TABV0072 

TABV0073 

TABV0074 

TABV0075 

TABV0076 

TABV0077 

TABV0078 

TABV0079 

TABV0080 

TABV0081 

TABV0082 

TABV0083 

TABV0084 

TABV0085 

TABV0086 

TABV0087 

TABV0088 

TABV0089 

TABV0090 

TABV0091 

TABV0092 

TABV0093 
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C  CHECK  FOR  TABLET  INITALIZED 

30  IFCINIT.NE.DGO  TO  40 
ICHAR=87 
XV=XB 
YVrYB 
INITrO 
RETURN 

40  CALL  BELL 

CALL  ONEPNT(IX,IY) 

C  CHECK  TO  SEE  IF  POINT  SENT  IS  A  MENU  COMMAND 

IF(IX.GT.MX2.0R.IX.LT.MX1)G0  TO  50 
IF(IY.GT.MY1.0R.n.LT.MY2)G0  TO  50 
IC=(IX-MX1)/LS+1 
IR=(MY1-IY)/LS+1 

C  CONVERT  RCW  AND  COLUMN  POSITION  TO  COMMAND  CHARACTER 

ICHAR=ICONV(IR,IC) 

IF(ICHAR.EQ.32)  return 

MSG4(18)=ICHAR 

LCNT=LCNT-LDEL 

C  LAST  MESSAGE  COMMAND 

CALL  NOTATE (0,LCNT,1 8, MSG4) 

IF(ICHAR.EQ.105)G0  TO  10 
IF ( IRETN ( IR , IC) .EQ. 1 ) RETURN 
GO  TO  40 

C  CONVERT  TABLET  UNITS  TO  VIRTUAL  UNITS  WITH  ANGLE  CORRECTION 

50  DX=IX-MXB 
DY=IY-MYB 

IF(DX.EQ.0,)DX=1.E-20 

R=SQRT(DX»DX+DY*DY) 

ANGR  =ATAN2 ( DY , DX ) -ANG 

XV=R«FACX*COS(ANGR )+XB 

YV=R*FACY*SIN(ANGR)+YB 

RETURN 

END 

CDRAWIT 

SUBROUTINE  DRAWIT ( NL , NPTA , X , Y , ISUB) 

COMMON/TKTRNX/ITEKC( 60 ) 

COMMON/TEKGPPR/DUM( 20 ) , EN ( 2 ) , DEL ( 2 ) , BEG ( 2 ) , RDX2 , RDY2 ,  ^LINE , 

1  NDRAW(30),MODE(30) 

DIMENSION  QSY( 306 ) , QSX( 306 ) , NPTA( 1) , X( 1 ) , Y( 1 ) , ISUB( 1 ) 
EQUIVALENCE  ( lOFF , ITEKC ( 30 ) ) 


0  SYMBOLS  1  LINE  2  SPLINE  WRT  X  3  SPLINE  WRT  Y  4  ARC  FIT  5  CLOSED 


TABV0094 

TABV0095 

TABV0096 

TABV0097 

TABV0098 

TABV0099 

TABV0100 

TABV0101 

TABV0102 

TABV0103 

TABV0104 

TABV0105 

TABV0106 

TABV0107 

TABV0108 

TABV0109 

TABV0110 

TABV0111 

TABV0112 

TABV0113 

TABV0114 

TABV0115 

TABV0116 

TABV0117 

TABV0118 

TABV0119 

TABV0120 

TABV0121 

TABV0122 

TABV0123 

TABV0124 

TABV0125 

TABV0126 

TABV0127 

DRAW0001 

DRAW0002 

DRAW0003 

DRAW0004 

DRAW0005 

DRAW0006 

DRAW0007 

DRAW0009 

DRAW0010 


IT=64 

K=1 

NSUM=1 


DRAW0012 

DRAW0013 

DRAW0014 

DRAW0015 
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NC=0 

DO  290  1=1,  NL 
ISY1^=MCDE(I) 

ITyP=IABS(ISYf'1) 

IF(ITYP.GT.I)  GO  TO  40 
NEND=NSUM+NPTA(I)-1 
IT=IT+1 

IF(IT.GT.90)IT=65 
DO  20  J=NSUM,NEND 
IF(NLINE.EQ.O)  GO  TO  10 
IF(NDR.'«V(I).EQ.O)  GO  TO  20 
10  XP=X(K) 

YP=Y(K) 

IF(J.EQ.NSUM)  CALL  MOVEA(XP,YP) 
IF(ITYP.EQ.I)  CALL  DRAWA(XP,YP) 
IF(ISYM.LT.O)  GO  TO  20 
CALL  MOVEA(XP,YP) 

IF(IOFF.EQ.O)  CALL  ANCHO(IT) 

CALL  MOVEA(XP,YP) 

20  K=ISUB(K) 

30  NSUM=NEND+1 
GO  TO  290 

PLOT  WITH  SPLINE 

40  NS=NC 

NPT=NPTA(I) 

NC=NC+NPT 
IT=IT+1 

IF(IT.GT.90)IT=65 
IF(NLINE,EQ.O)  GO  TO  60 
IF(NDRAW(I).NE.O)  GO  TO  60 
C  LOCATE  POINTER  OT  NEXT  LINE 

DO  50  L=1,NPT 
50  K=ISUB(K) 

GO  TO  290 
60  JFIT=2 
YO=Y(K) 

K1=ISUB(K) 

IF(ITYP.GT.2)  GO  TO  80 
XO=X(K) 

C  CHECK  X  DATA  FOR  ASCENDING  ORDER 

DO  70  L=2,NPT 
X1=X(K1 ) 

IF(XI.LE.XO)  GO  TO  110 
K1=ISUB(K1) 

70  XO=X1 
GO  TO  210 


DRAW0016 


DRAW0008 

DRAW0017 

DRAW0018 

DRAW0019 

DRAW0020 

DRAW0021 

DRAW0022 

DRAW0023 

DRAW0024 

DRAW0025 

DRAW0026 

DRAW0027 

DRAW0028 

DRAW0029 

DRAWOO 

DRAW0030 

DRAW0031 

DRAW0033 

DRAW0034 

DRAW0035 

DRAW0036 

DRAW0042 

DRAW0043 

DRAW0044 

DRAW0045 

DRAW0046 

DRAW0047 

DRAW0048 

DRAW0049 

DRAW0050 

DRAW0051 

DRAW0052 

DRAW0053 

DRAW0054 

DRAW0055 

DRAW0056 

DRAW0057 

DRAW0058 

DRAW0059 

DRAW0060 

DRAW0061 

DRAW0062 

DRAW0063 
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80  IF(ITYP.GT.3)  GO  TO  100 

DRAWC055 

CHECK  Y  DATA  FOR  ASCENDING  ORDER 

DRAWCO60 

DO  90  L=2,NPT 

DRAWC067 

Y1=Y(K1) 

DRAWCC68 

IF(YI.LE.YO)  GO  TO  110 

DRAW0069 

K1=ISUB(K1) 

DRAW0070 

90  YO=Y1 

DRAWG071 

GO  TO  210 

DRAW0072 

100  JFIT=ITYP-2 

DRAW0073 

110  NCIR=0 

DRAWOO74 

IF(JFIT.EQ.3)  NCIR=-NPT/2-1 

DRAWQ075 

MPT=NPT-2*NCIR 

DRAW0076 

QSY(1)=MPT 

DRAWC077 

QSX(1)=MPT 

DRAW0078 

S=0. 

DRAW0079 

KA=NS 

DRAWOOSO 

KO=KA 

DRAWC081 

KE=KO+NPT 

DRAW0082 

KSAVErK 

DRAW0083 

KA=KA+NCIR 

DRAW0084 

DO  160  M=1,MPT 

DRAW0085 

M1=M+1 

DRAW0086 

KA=KA+1 

DRAW0087 

IF(KA.GT.KO)  GO  TO  130 

DRAW0088 

NDOsNPT+NCIR 

DRAW0089 

DO  120  11=1, NDO 

DRAWOO9O 

120  K=ISUB(K) 

DRAWOO9I 

KA=KA+NPT 

DRAWCO92 

GO  TO  140 

DRAWOO93 

130  IF(KA.NE.(KE+1))GO  TO  140 

DRAWOO94 

JSAVE=K 

DRAWOO95 

KrKSAVE 

DRAW0096 

KA=KA-NPT 

DRAWOO97 

140  CONTINUE 

DRAW0098 

L=M1+MPT 

DRAWOO99 

YYYP=Y(K) 

DRAW0100 

XXXP=X(K) 

DRAW0101 

K=ISUB(K) 

DRAW0102 

IFCM.EQ. 1)  GO  TO  150 

DRAWOIO3 

DS=SQRT ( RDX2* ( XXXP-XO ) **2+RDY2* ( YYYP-YO ) *»2 ) 

DRAW0104 

S=S+DS 

DRAW0105 

150  XO=XXXP 

DRAW0106 

YOrYYYP 

DRAW0107 

QSX(M1)=S 

DRAW0108 

QSY(M1)=S 

DRAWOIO9 

QSX(L)=XXXP 

DRAW0110 

160  QSY(L)=YYYP 

DRAWOni 

KA=K0+NPT 

DRAW0112 
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QSX(L-(-1)=0.  DRAW0113 
QSY(U1)=0.  DRAW0114 
QSX(L-t-2)=1.  DRAW0115 
QSY(L-i-2)=1.  DRAW0116 
X0=QSX(MFT-i-2-NCIR )  DRAW01 1 7 
Y0=QSY(MFT-(.2-NCIR )  DRAW01 1 8 
CALL  MOVEA(XO,TO)  DRAW0119 
IF(IOFF.EQ.O)  CALL  ANCHO(IT)  DRAW0120 
CALL  MOVEA(XO,YO)  DRAW0121 
SCK=QSX(3-NCIR)  DRAW0122 
S=QSX(2-NCIR)  DRAW0123 
IF(NPT.LE.I)  GO  TO  290  DRAW0124 
DC =40.  DRAW0125 
DS=40.  DRAy0126 


NCK=2 
170  S=S-hDS 

XP=SPLNQ1(1,QSX,S) 

YP=SPLNQ1(1,QSY,S) 

DCK=SQRT(RDX2»(XO-XP)**2+RDY2*(YO-YP)**2) 

DS=  DC*DS/DCK 
180  IF(S.LT.SCK)  GO  TO  200 
NSYM=MFTH-1-t-NCK-NCIR 
XS=QSX(NSYM) 

YS=QSY(NSYM) 

CALL  DRAWA(XS,YS) 

IF(ISYM.LE.O.AND.NCK.NE.NPT)  GO  TO  190 
CALL  MOVEA(XS,YS) 

IF(IOFF.EQ.O)  CALL  ANCHO(IT) 

CALL  MOVEA(XS,YS) 

190  NCK=NCK-i-1 

SCK=QSX(NCK-f1-NCIR) 

IF(NCK.LE.NPT-t-JFIT-2)  GO  TO  180 
IF(JFIT.EQ.3)  K=JSAVE 
GO  TO  290 

200  CALL  DRAWA(XP,YP) 

XO=XP 
YO=YP 
GO  TO  170 
210  QSX(1)=NPT 

DO  240  M=1,NPT 

N=M-k1 

KA=NS+M 

L=N-hNPT 

XP=X(K) 

YP=Y(K) 

IF(M.NE.I.AND.M.NE.NPT.AND.ISYM.LE.O)  GO  TO  220 
CALL  MOVEA(XP,YP) 

IF(IOFF.EQ.O)  CALL  ANCHO(IT) 


DRAW0127 
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DRAW0129 

DRAW0130 
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DRAW0149 

DRAW0150 

DRAW0151 

DRAW0152 

DRAW0153 
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220  IFCITYP.rJE.S)  GC  TO  230 

DRAW0161 

QSX(N)=YP 

DRAW0162 

QSX(L)=XP 

DRAW0163 

GO  TO  240 

DRAW0164 

230  QSX(N)=XP 

DRAW0165 

QSX(L)=YP 

DRAW0166 

240  K=ISUB(K) 

DRAW0167 

QSX(L+1)=0. 

DRAW0168 

QSX(L+2)=1. 

DRAW0169 

XEN=QSX(NPT+1 ) 

DRAW0170 

XIN=QSX(2) 

DRAW0171 

IFITP=ITYP-1 

DRAW0172 

BCK=BEG(IFITP) 

DRAW0173 

ECK=EN(IFITP) 

DRAW0174 

DELT=DEL(IFITP)*30. 

DRAW0175 

IF(XIN.LT.BCK)  XIN=BCK 

DRAW02176 

IF(XEN.GT.ECK)  XEN=ECK 

DRAW0177 

KILLrO 

DRAW0178 

DO  280  Mr 1,200 

DRAW0179 

XIrXIN+DELT*(M-1 ) 

DRAW0180 

IF(XI.LT.XEN)  GO  TO  250 

DRAW0181 

KILLrl 

DRAW0182 

XIrXEN 

DRAW0183 

250  YIrSFLNQ1(1,QSX,XI) 

DRAW0184 

IF(ITYP.EQ.3)  GO  TO  260 

DRAW0185 

XPrXI 

DRAW0186 

YPrYI 

DRAW0187 

GO  TO  270 

DRAW0188 

260  XPrYI 

DRAW0189 

YPrXI 

DRAW0190 

270  IF(M.EQ.I)  CALL  MOVEA(XP,YP) 

DRAW0191 

CALL  DRAWA(XP,YP) 

DRAW0192 

IFCNPT.EQ. 1)  GO  TO  290 

DRAW0193 

IF(KILL.EQ.I)  GO  TO  290 

DRAW0194 

280  CONTINUE 

DRAW0195 

290  CONTINUE 

DRAW0196 

300  RETURN 

DRAW0197 

END 

DRAW0198 

CSPLNQ1 

FUNCTION  SPLNQ1  (NLOC ,X,XINDEP) 

SPLN0001 

C***  LOCAL  CUBIC  FIT  8/9/77  M.J.  CADDY 

SPLN0002 

DIMENSION  X(1),QM(3) 

SPLN0003 

EQUIVALENCE  (QM( 1) ,T3) , (QM(2) ,Q2) , (QM(3) ,Q3) 

SPLN0004 

XINrXINDEP 

SPLN0005 

NSrNLOC 

SPLN0006 

NOPTSrX(NS) 

SPLN0007 

IDrNS+NOPTS 

SPLN0008 

NSP1rNS+1 

SPLN0009 
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NSP2=NS+2 

SPLN0010 

IF(NOPTS.LE.I)  GO  TO  130 

SPLN001 1 

IF(N0PTS.GT.2)  GO  TO  10 

SPLN0012 

M=ID+NOPTS 

SPLN0013 

T3=(X(N)-X(N-1 ))/(X(ID)-X(ID-1)) 

SPLN0014 

MrlD 

SPLN0015 

NTRAP=1 

SPLN0016 

GO  TO  280 

SPLN0017 

10  NS2=N0PTS*2+NSP1 

SPLN0018 

L=X(NS2) 

SPLN0019 

LSC=NS2+1 

SPLN0020 

IQMODE=X(LSC) 

SPLN0021 

K=L+NS 

SPLN0022 

NLrNSPI 

SPLN0023 

NHrlD 

.  SPLN0024 

NTRAP=-1 

SPLN0025 

C***  binary  search  for  interval 

SPLN0026 

IF(XIN-X(ID))30, 140,20 

SPLN0027 

20  NTRAPrO 

SPLN0028 

GO  TO  150 

SPLN0029 

30  IF(XIN-X(NSP1))40,40,60 

SPLN0030 

40  NTRAPrI 

SPLN0031 

50  K=NSP2 

SPLN0032 

GO  TO  160 

SPLN0033 

60  IF.(L)  120, 120,70 

SPLN0034 

70  IF(XIN-X(K))80,100,100 

SPLN0035 

80  NH=K 

SPLN0036 

K=K-1 

SPLN0037 

90  IF(XIN-X(K))110,100,100 

SPLN0038 

100  NL=K 

SPLM0039 

GO  TO  120 

SPLN0040 

110  NH=K 

SPLN0041 

120  K=(NH-NL)/2+NL 

SF«N0042 

IF(K-NL)90, 140,90 

SPLN0043 

130  YOUT=X(NSP2) 

SPLN0044 

GO  TO  320 

SPLN0045 

140  LFAST=L-NH+NS 

SPLN0046 

X(NS2)=NH-NS 

SPLN0047 

150  K=NH 

SPLN0048 

160  M=K 

SPLN0049 

N=M+NOPTS 

SPLN0050 

Y3=X(N-1) 

SPLN0051 

X3=X(M-1) 

SPLN0052 

C***  CHECK  FOR  FAST  MODE  AND  EXTRAPOUTION 

SPLN0053 

IF(NTRAP.GE.O)  GO  TO  180 

SPLN0054 

IF(IQMODE*L.EQ.O.OR.LFAST.NE.O)  GO  TO  180 

SPLN0055 

DO  170  1=1,3 

SPLN0056 

170  QM(I)=X(LSC+I) 

SPLN0057 
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V 


GO  TO  3T0 

5FL:]C053 

180  Y4=X(N) 

spl:ico59 

X4=X(M) 

A3=X4-X8 

3Fl:ico6i 

S3=(Y4-Y3)/A3 

5PLMCC62 

IF(M.EQ.NSP2)  GO  TO  190 

5PLNC063 

X2=X(M-2) 

3FLN0C64 

Y2=X(N-2) 

3PLN0065 

S2=(Y3-Y2)/(X3-X2) 

3PLN0066 

IF(M.EQ.ID)  GO  TO  200 

5PLN0067 

190  X5=X(M-h1) 

5PLN0068 

Y5=X(N-h1) 

3PLN0069 

S4=(Y5-Y4)/(X5-X4) 

5R.N0070 

IF(M.EQ.NSP2)  S2=S3-f53-S4 

5R-N0071 

GO  TO  210 

5PLN0072 

200  S4=S3-i-S3-S2 

3F-N0073 

210  IF(M.LE.(NSP2-h1))  GO  TO  220 

5PLN0074 

S1=(Y2-X(N-3))/(X2-X(M-3)) 

5PLN0075 

GO  TO  230 

3PLN0076 

220  51=52-1-32-53 

3PLN0077 

230  IF(M.GE.(ID-1))  GO  TO  240 

5PLN0078 

35=  (X  (N-k2  ) -Y5 )  /  ( X  (M-h2  ) -X5 ) 

3R.N0079 

GO  TO  250 

5F.N0080 

240  55=54+34-33 

3PLN0081 

250  W2=AB5 (34-33) 

3PLN0082 

W3=AB5 (52-31) 

3PLN0083 

3W=W2+W3 

3PLN00S4 

IF(3W.NE.0.0)  GO  TO  260 

5PLN0085 

W2=0.5 

3PLN0086 

W3=0.5 

5PLN0087 

3W=1.0 

3n-N0088 

260  T3=(W2*32+W3*33)/5W 

5PLN0089 

W3=AB3 (55-34) 

5PLN0090 

W4=AB3 (33-32) 

3PLN0091 

5W=W3+W4 

•3PLN0092 

IF(3W.NE.0.0)  GO  TO  270 

SPLN0093 

W3=0.5 

3PLN0094 

W4=0.5 

3PLM0095 

5W=1.0 

5PLN0096 

270  T4=(W3*53+W4*34)/5W 

3F-N0097 

IF(NTRAP.LT.O)  GO  TO  290 

3PLN0098 

IF(NTRAP.EQ.O)  T3=T4 

SPLN0099 

280  IX=M-NTRAP 

3PLN0100 

C»»*  FA3T  EXIT  FOR  2  POINT3  AND  LINEAR  EXTRAPOLATION 

5PLN0101 

YOUTrX ( IX+NOPT3 ) + (XIN-X ( IX) ) »T3 

3PLN0102 

GO  TO  320 

SPLN0103 

290  Q2=(2.0*(S3-T3)+33-T4)/A3 

3PLN0104 

Q3=(-S3-S3+T3-^T4)/(A3*A3) 

SPLN0105 
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IF(IQMODE»LFAST.EQ.O)  (30  TO  310 
DO  300  1=1,3 
300  X(LSC+I)=QM(I) 

310  DX=XIN-X3 

YOUT=Y3+DX»(T3+DX*(Q2+DX»Q3 ) ) 

320  SPLNQIrYOUT 
RETURN 
END 

CTTITE 

SUBROUTINE  TTITE(IX,IY,NTL,LABTL,NM,IA) 
DIMENSION  LABTL(1),IP(136) 

NIL  rNUMBER  OF  10  CHARACTER  WORDS 
NM  MAX  CHARACTERS  PER  LINE 
lA  SWITCH, IA=0  H0RIZ,IA=1  VERTICAL 
IX  SCREEN  CENTER 
lY  SCREEN  CENTER 
IF(NTL.LE.O)  RETURN 
NC=10*NTL 

C  GET  CHARACTER  SIZE 

CALL  CSIZE(IHORZ,IVERT) 

C  CONVERT  LABEL  TO  ADE 

CALL  KAM2AS(NC,UBTL,IP) 

1X1 =IX 
IY1=IY 
ITL1=0 
NBLKrO 

DO  70  K=1,NC 

C  CHECK  FOR  LEADING  BLANKS 

IF(IP(K).NE.32)  GO  TO  10 
IF(ITLl.EQ.O)  GO  TO  70 
NBLK=NBLK+1 

C  CHECK  FOR  3  BLANKS  TO  TERMINATE  LINE 

IF(NBLK.NE.3)  GO  TO  20 
ITL1=ITL1-2 
GO  TO  50 
10  NBLKrO 

C  CHECK  FOR  MAX  LINE  LENGTH  EXCEEDED 

20  IF(ITLl.LT.NM)  GO  TO  30 
IF(IP(K).EQ.32)  GO  TO  50 
30  ITL1=ITL1+1 
IP(ITL1)=IP(K) 

IF(K.LT.NC)  GO  TO  70 
40  ITL1=ITLT-NBLK 

C  CHECK  FOR  VERTICAL  OR  HORIZ  LABEL 

50  IF(IA.NE.O)  GO  TO  60 
IX1=IX-IH0RZ*ITL1*.5 
CALL  NOTATE ( 1X1, IY1,ITL1, IP) 
IY1=IY1-IVERT*1.1 


SPLN0106 

SPLN0107 

SPLN0108 

SPLN0109 

SPLN0110 

SPLN0111 

SPLN0112 

SPLN0113 

TTIT0001 

TTIT0002  , 

TTIT0003 

TTIT0004 

TTIT0005  - 

TTIT0006 

TTIT0007 

TTIT0008 

TTIT0009 

TTIT0010 

TTIT0011 

TTIT0012 

TTIT0013 

TTIT0014 

TTIT0015 

TTIT0016 

TTIT0017 

TTIT0018 

TTIT0019 

TTIT0020 

TTIT0021 

TTIT0022 

TTIT0023 

TTIT0024 

TTIT0025 

TTIT0026 

TTIT0027 

TTIT0028 

TTIT0029 

TTIT0030  > 

TTIT0031 

TTIT0032 

TTIT0033  > 

TTIT0034 

TTIT0035 

TTIT0036 

TTIT0037 

TTIT0038 

TTIT0039 

TTIT0040 
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ITLIrO 
GO  TO  70 

60  IY1=IY+IVERT»ITL1*.5 
CALL  MOVABS(IXl,ni) 
CALL  VLABEL(ITL1,IP) 
IX1=IX1+IH0RZ*1.1 
ITLIrO 
70  CONTINUE 
END 

CGETIN 

SUBROUTINE  GETIN(NIN,Y) 


MICHAEL  CADDY  3/19/78 
DIMENSION  Y(1),IC(80) 

C  FREE  FORM  INPUT  CODE 

NWrO 
10  JrO 

READ  20, IC 

C  CHECK  FOR  END  OF  FILE 

IF(E0F(5).EQ.0)G0  TO  30 
NINrNW 
RETURN 

20  FORMAT(80R1) 

30  JCrO 
JD=0 

JSrI 

NCrO 

XrO. 

40  JrJ+1 

C  ONLY  ONE  CARD  PER  INPUT  READ 

C  MODIFIED  TO  READ  MORE  THAN  ONE  CARD  4/26/78  MJC 

IF(J.GT.80)  GO  TO  10 
IrlC(J) 

C  CHECK  FOR  VALID  NUMERIC  FIELD 

IF(I.GT.32B.AND.I.LT.45B)  GO  TO  110 
C  IGNORE  LEAD  +  SIGN 

IF(I.EQ.45B)  GO  TO  40 
C  SET  FLAG  FOR  NEGATIVE  VALUE 

IF(I.NE.46B)  GO  TO  50 
JSr-1 
GO  TO  40 

C  CHECK  FOR  DECIMAL 

50  IF(I.NE.57B)  GO  TO  60 
IF(JC.EQ.-I)  GO  TO  120 

C  IF  THIS  IS  SECOND  DECIMAL  BLOW  OFF  TO  ERROR  CODE 

JCr-1 
GO  TO  40 


TTIT0041 

TTIT0042 

TTIT0043 

TTIT0044 

TTIT0045 

TTIT0046 

TTIT0047 

TTIT0048 

TTIT0049 

GETI0001 

GETI0002 

GETI0003 

GETI0004 

GETI0005 

GETI0006 

GETI0007 

GETI0008 

GETI0009 

GETI0010 

GETI0011 

GETI0012 

GETI0013 

GETI0014 

GETI0015 

GETI0016 

GETI0017 

GETI0018 

GETI0019 

GETI0020 

GETI0021 

GETI0022 

GETI0023 

GETI0024 

GETI0025 

GETI0026 

GETI0027 

GETI0028 

GETI0029 

GETI0030 

GETI0031 

GETI0032 

GETI0033 

GETI0034 

GETI0035 

GETI0036 

GETI0037 

GETI0038 
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C  CHARACTER  IS  BLANK  TREAT  AS  COMMA  IF  NOT  LEADING 

60  IF(I.EQ.55B)GO  TO  80 
IF(I.EQ.56B)GO  TO  70 
GO  TO  120 

70  IF(NC.GT.O)GO  TO  90 
NWrNW-t-l 

IF(NW.Gr.NIN)  RETURN 
GO  TO  HO 

C  TWO  COMMAS  ..IGNORE  THIS  DATA  FIELD  AND  GO  ON  TO  NEXT 

80  IF(NC.EQ.O)  GO  TO  40 
C  SHIFT  DECIMAL  TO  NUMBER 

90  X=JS*X*10.**JD 
NW=NW-k1 
Y(NW)=X 

if(nw.ge.nin)  return 

GO  TO  30 
110  JDrJD-t-JC 
NCrNC-t-l 

C  ADD  DIGIT  TO  NUMBER  ,, CAREFULLY 

X=X»10-h(I-33B) 

GO  TO  40 

C  ERROR  CODE 

120  DO  130  K=1,80 

130  IC(K)=55B 
IC(J)=47B 
PRINT  140  ,IC 
140  F0RMAT(2X,80R1) 

PRINT  150 

150  FORMAT(*  BAD  FIELD,  RE-ENTER  DATA*) 

GO  TO  10 
END 


GETI0039 

GETI0040 

GETI0041 

GETI0042 

GETI0043 

GETI0044 

GETI0045 

GETI0046 

GETI0047 

GETI0048 

GETI0049 

GETI0050  , 

GETI0051 

GETI0052 

GETI0053 

GETI0054 

GETI0055 

GETI0056 

GETI0057 

GETI0058 

GETI0059 

GETI0060 

GETI0061 

GETI0062 

GETI0063 

GETI0064 

GETI0065 

GETI0066 

GETI0067 

GETI0068 

GETI0069 
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